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PREFACE 

This  report  describes  a  computerized  method  for  analysis  and 
design  of  pile  groups  that  is  currently  being  used  by  several  Corps  of 
Engineers  offices.  Criteria  for  a  new,  comprehensive  computer  program 
for  pile  analysis  and  design  are  also  discussed.  The  work  was  sponsored 
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CONVERSION  FACTORS,  INCH-FOUND  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


Inch- pound  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 


B)L 


To  Obtain 


feet 

inches 

inch-pounds  (force) 
pounds  (force)  per  inch 
pounds  (force)  per  square  foot 
pounds  (force)  per  square  inch 
pounds  (mass)  per  cubic  foot 
pounds  (mass)  per  cubic  inch 


0.3048 

2.54 

0.1129848 

0.1751268 

47.880263 

6.894757 

16.01846 

0.0276799 


metres 

centimetres 

newton-metres 

kilonewtons  per  metre 

pascals 

kilopascals 

kilograms  per  cubic  metre 

kilograms  per  cubic 
centimetre 


tons  (2000  lb  mass) 


907.18474  kilograms 
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SCOPE 


The  purpose  of  this  paper  is  to  present  one  method  for  pile  group 
design  and  analysis  as  practiced  by  the  Corps  of  Engineers,  and  to  propose 
criteria  for  systematizing  this  method  in  a  computer  program.  This  paper 
describes  a  computerized  method  of  pile  group  analysis  (including  sample 
problems)  and  lists  criteria  for  a  new,  more  comprehensive  program;  it 
includes  an  overview  of  advanced  methods  of  pile  design,  and  briefly 
discusses  selection  of  pile  types,  methods  of  installation,  and  allowable 
stresses 

II.  PILE  DESIGN  CONSIDERATIONS. 

A.  Economic .  Pile  foundations  are  a  major  cost  in  a  structure.  Pile 
foundations  that  provide  the  lowest  first  cost  are  of  paramount  importance. 

A  cost  comparison  must  be  made  of  the  relative  cost  of  different  type  piles 
and  cost  of  installation.  Scheduling  and  availability  may  affect  pile 
costs.  Details  affecting  selection  of  pile  type  are  presented  in  Appendix  A. 

B.  Affect  on  Adjacent  Structures.  Proximity  of  adjacent  structures 
may  dictate  the  type  of  pile  or  installation  used.  Adverse  effects  of  soil 
displacement  or  vibration  caused  by  driving  piles  may  compel  the  use  of 
drilled  caissons,  nondisplacement  piles,  or  jetting  or  predrilling  of  piles. 

C.  Difficulty  in  Installation.  Hard  strata,  boulders,  buried  debris, 
and  other  obstructions  may  necessitate  the  use  of  piles  durable  enough  to 
sustain  driving  stresses.  Jetting,  predrilling,  or  spudding  may  be 
required.  Descriptions  of  installation  methods  are  presented  in  Appendix  B. 

D.  Environment.  Corrosion  in  sea  water  will  require  consideration  for 
protective  coating,  concrete  jacketing,  or  cathodic  protection  if  steel 
piling  is  used.  The  presence  of  marine-borers  may  negate  the  use  of  wood 
piling  and  the  subsequent  use  of  steel  or  concrete  piling. 


S.  Displacements.  Limitations  on  lateral  or  rotational  movement  will 
affect  tbe  type  of  pile  used  and  the  configuration  of  the  pile  group. 

Stiffer  piles  and  the  degree  of  fixity  to  the  pile  caps  are  considerations  to 
limit  displacements. 

F.  Foundation  Materials.  The  capacity  of  the  piling  may  be  limited  by 
failure  of  the  foundation  materials,  evidenced  by  excessive  settlement  of 
piles  under  applied  load.  The  capacity  of  the  foundation  materials  is 
usually  evaluated  by  static  resistance  formulas  during  design  and  verified  by 
load  tests  prior  to  construction.  Dynamic  driving  formulas  are  generally  not 
a  reliable  basis  for  estimating  pile  capacities  unless  correlated  with  load 
tests  and  previous  experience  at  similar,  nearby  sites.  More  reliable 
predictions  of  dynamic  behavior  during  driving  are  based  on  complex 
computerized  models  of  hammer-pile-soil-  interaction  using  the  ID  wave 
equation. 

G.  Failure  Modes. 

1.  Bearing  capacity  failure  of  the  pile-soil  system. 

2.  Excessive  settlement  due  to  compression  and  consolidation  of 
the  underlying  soil. 

3.  Structural  failure  of  the  pile  under  service  loads. 

4.  Bearing  capacity  failure  caused  by  improper  installation 

methods . 

5.  Structural  failure  resulting  from  detrimental  pile 
installation.  This  may  be  due  to  unforseen  subsoil  conditions  or  to  freeze, 
compaction,  liquification,  or  heave  of  the  soil.  It  could  be  caused  by 
driving  sequence,  size  of  hammer,  vibration,  over  or  under  driving,  improper 
preexcavation  methods,  substitution  of  materials,  improper  workmanship,  or 
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limitations  of  the  Contractor*}  equipment  or  expertise.  These  conditions  are 
described  in  detail  in  Chapter  2  of  Reference  1. 

H.  Other  Considerations.  For  more  detailed  discussion  of  the  above 
considerations,  and  for  others  not  mentioned  above,  see  Reference  2. 

III.  BASIC  PILE  GROUP  ANALYSIS. 

This  section  presents  the  fundamentals  of  a  basic  method  of  pile  group 
analysis  which  is  currently  available  in  various  computer  programs  including 
LMVDPILE  which  is  in  WESLIB.  Several  hand  analysis  methods  are  shown  with 
the  sample  problems  in  Appendix  E.  This  computer  method  is  capable  of 
handling  three-dimensional  loading  and  pile  geometry.  It  is  valid  for  static 
analysis  of  a  linear,  elastic  system.  Interaction  between  pile  and  structure 
is  limited  to  the  extremes  of  a  fully  fixed  or  fully  pinned  connection. 
Interaction  between  the  pile  and  soil  is  represented  by  a  linear,  elastic 
pile  stiffness  (applied  load  per  unit  deflection)  at  the  top  of  the  pile. 

The  base  of  the  structure  is  assumed  to  act  as  a  rigid  body  pile  cap 
connecting  all  piles;  the  cap  flexibility  is  not  considered. 

A.  Basic  Analysis  Method.  The  basic  pile  group  analysis  method 
represents  each  pile  by  its  calculated  stiffness  coefficient,  in  the  manner 
proposed  by  laul  (3).  The  stiffness  coefficients  of  all  piles  are  summed  to 
determine  a  stiffness  matrix  for  the  total  pile  group.  Displacements  of  the 
rigid  pile  cap  are  determined  by  multiplying  the  sets  of  applied  loads  by  the 
inverse  of  the  group  stiffness  matrix.  Displacements  of  the  rigid  pile  cap 
define  deflections  of  individual  pile  heads  which  are  then  multiplied  by  the 
pile  stiffness  coefficients  to  determine  the  forces  acting  on  each  pile 
head.  The  key  step  in  the  method  is  in  determining  individual  pile  stiffness 
coefficients,  at  the  pile  head,  based  on  known  or  assumed  properties  of  pile 
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and  soil.  Since  this  is  a  three-dimensional  analysis  method,  each  pile  head 
has  six  degrees  of  freedom  (DOF),  three  translations  and  three  rotations.  A 
stiffness  coefficient  must  be  determined  for  each  DOF  and  for  all  coupling 
effects  (e.g.  lateral  deflection  due  to  applied  moment).  The  pile  location 
and  batter  angle  are  also  accounted  for  when  individual  pile  stiffness 
coefficients  are  combined  to  form  the  total  stiffness  matrix  for  the  pile 
group. 

B.  Pile-Structure  Interaction.  Piles  are  mathematically  represented 
in  the  analysis  by  their  axial,  lateral  and  rotational  stiffness,  as  springs 
resisting  motion  of  the  rigid  cap.  Such  a  system  is  shown  in  Figure  1. 


As  mentioned  above,  consideration  is  given  only  to  piles  which  are  fully 
fixed  or  pinned  to  the  pile  cap.  A  pile  embedded  only  a  short  distance  into 
the  cap  may  be  assumed  to  transfer  no  moment  at  the  pile  head.  Such  a  pile 
will  resist  only  shear  and  axial  loads.  Well-embedded  piles  will  resist 
shears,  moments,  and  axial  loads  and  will  have  coupling  stiffness,  referred 
to  above.  It  is  necessary  to  consider  the  fixity  of  the  cap-pile  joint  to 
adequately  determine  pile  stiffness.  This  should  be  done  in  conjunction  with 
consideration  of  pile-soil  interaction.  Pile  head  fixity  parameters  have 
been  derived  by  Dawkins  (4).  Once  an  analysis  has  determined  the  forces 
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acting  on  each  pile,  theae  forces  may  then  be  applied  to  the  pile  cap  to 
determine  its  internal  shears  and  moments.  However,  analysis  of  the  pile  cap 
is  outside  the  scope  of  this  section. 

C.  Pile-Soil  Interaction.  Interaction  between  the  pile  and  soil  is 
the  most  important  consideration  in  determining  pile  stiffness.  Therefore, 
it  is  necessary  to  have  reliable  information  about  soil  properties.  Soil 
properties  can  affect  the  axial,  lateral,  or  torsional  stiffness  of  the 
pile.  The  type  of  loading  expected  (static  or  cyclic)  and  the  pile  spacing 
should  also  be  considered  since  cyclic  loading  or  close  spacing  may  both 
reduce  individual  pile  stiffness. 

1.  Axial  Stiffness.  Axial  load  in  a  pile  may  be  transferred  to 
the  soil  by  some  combination  of  tip  bearing  and  skin  friction.  For  a  pile 
transferring  all  load  by  tip  bearing  the  axial  stiffness  is  obviously  AE/L, 
the  stiffness  of  any  axially  loaded  structural  member.  For  a  pile 
transferring  all  load  by  skin  friction  uniformly  along  its  entire  length, 
with  no  further  tip  movement,  the  axial  stiffness  is  2AE/L.  Any  other  axial 
stiffness  is  possible  for  other  types  of  piles,  as  shown  in  Figure  2. 


Axial 

Stiffness 

Coefficient: 


T 


Tip  Bearing 
P 


AE/L 


f 

i  \ 
i  \ 
\  f 


No  Penetration  f  Tip  Penetration 


2  AE/L  K  AE/L 

FIGURE  2 
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A  further  complication  of  pile  axial  stiffness  involves  consideration  of 
tension  piles.  Generally,  a  pile  in  tension  will  be  less  stiff  than  the  same 
pile  in  compression.  Since  only  a  single  elastic  stiffness  coefficient  may 
be  specified  for  each  pile,  that  stiffness  must  be  based  on  whether  the  load 
is  expected  to  be  tension  or  compression. 

2.  Lateral  Stiffness.  Pile  lateral  stiffness  refers  to 
rotational  stiffness  and  coupling  effects,  in  addition  to  actual 
translational  stiffness.  The  most  important  consideration  is  the  resistance 
of  the  soil  to  translation  of  a  pile.  The  degree  of  fixity  between  the  cap 
and  the  pile  must  also  be  considered.  The  pile  may  be  represented  as  a  beam 
on  elastic  foundation,  with  the  soil  represented  as  a  set  of  springs  acting 
on  the  pile,  as  shown  in  Figure  3. 


Load 


•H Mr 
,hAAr 
H/W" 


— \A/V“£ 
■VW^ 


Soil  Stiffness 


FIGURE  3 


Though  soil  properties  are  often  highly  non-linear,  an  approximate  linear 
lateral  stiffness  coefficient  must  be  determined.  Several  analytical  methods 
may  be  used  to  determine  this  stiffness.  One  method  is  to  use  any  beam 
analysis  computer  program  capable  of  representing  the  beam-spring  system 
shown  in  Figure  3.  The  stiffness  equals  the  force  required  to  cause  a  unit 
displacement  at  the  pile  head.  This  method  may  be  used  to  determine  lateral, 
rotational  and  coupling  stiffness  coefficients.  A  method  for  determining 
appropriate  values  for  the  stiffness  of  the  soil  springs  is  included  as 
Appendix  D.  Methods  for  determining  pile  stiffness  coefficients  are 
presented  in  detail  in  Appendix  C. 

3.  Torsional  Stiffness.  For  groups  of  piles,  torsion  on 
individual  piles  is  usually  unimportant  and  may  be  neglected  by  using  zero 
torsional  stiffness.  Where  torsion  of  individual  piles  is  important, 
torsional  stiffness  may  be  determined  in  a  manner  similar  to  that  described 
above  for  axial  stiffness  (5). 

D.  Analysis  Details.  As  mentioned  above,  this  analysis  method  has 
been  systematized  for  use  in  computer  programs.  Several  of  these  programs 
were  identified  during  the  Corps-Wide  Conference  on  Computer  Aided  Design  in 
Structural  Engineering  (6).  Following  are  some  of  the  detailed  formulations 
used  in  these  programs. 

1.  Coordinate  System.  The  basic  coordinate  system  is  a  right 
hand  system,  as  shown  in  Figure  4.  The  three  axes  are  labeled  1,  2,  and  3, 
with  the  3  axis  being  positive  downward.  This  global  coordinate  system  is 
used  for  specification  of  pile  locations  and  orientations,  applied  forces  and 
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moments  on  the  pile  cap,  and  for  calculation  of  total  pile  group  stiffness 
and  resulting  pile  cap  displacements. 


Global  Coordinate  System 
FIGURE  4 


Each  pile  also  has  its  own  local  coordinate  system  as  shown  in  Figure  5.  The 
axes  are  labeled  1,  2,  and  3  and  are  located  by  specifying  translations  and 
rotations  from  the  global  coordinate  system.  The  3  axis  is  positive  along 
the  pile  length,  the  1  and  2  axes  correspond  to  the  pile  principal  axes,  and 
the  pile  batter  is  in  the  local  1-3  plane. 


'iP 

Global 

Coordinate 

System 


Local  Coordinate  System 


FIGURE  5 
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The  local  coordinate  system  is  used  for  calculation  of  the  stiffness 
coefficients,  displacements,  and  forces  of  individual  piles. 

2.  Pile  Stiffness  Matrix.  For  a  pile  with  6  degrees  of  freedom, 
individual  pile  stiffness  coefficients  are  represented  by  a  6  x  6  matrix: 


Where  1,  2,  and  3  refer  to  the  pile  coordinate  system  axes,  and  4,  5,  and  6 
are  rotations  about  those  axes.  Thus,  b^  and  b ^  are  lateral 
stiffnesses,  b^  is  axial  stiffness,  b^  and  b^  ere  rotational 
stiffnesses,  b66  is  torsional  stiffness,  and  b15,  b^,  b^2,  and  b^ 
are  coupling  stiffnesses. 

bjj  -  is  the  force  required  to  displace  the  pile  head  a  unit  distance 
along  the  local  1  axis 

^22  ~  *s  t*ie  f°rce  required  to  displace  the  pile  head  a  unit  distance 

along  the  local  2  axis 

b^  -  is  the  force  required  to  displace  the  pile  head  a  unit  distance 
along  the  local  3  axis 

b/(/)  -  is  the  moment  required  to  displace  the  pile  head  a  unit  rotation 

around  the  local  1  axis 
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55 


*b 


15 


*b 


24 


*b 


51 


*b 


42 


is  the  moment  required  to  displace  the  pile  head  a  unit  rotation 
around  the  local  2  axis 

is  the  force  along  the  local  1  axis  caused  by  a  unit  rotation  of 
the  pile  head  around  the  local  2  axis 

is  the  force  along  the  local  2  axis  caused  by  a  unit  rotation  of 
the  pile  head  around  the  local  1  axis 

is  the  moment  around  the  local  2  axis  caused  by  a  unit  displacement 
of  the  pile  head  along  the  local  1  axis 

is  the  moment  around  the  local  1  axis  caused  by  a  unit  displacement 
of  the  pile  head  along  the  local  2  axis 


*Since  the  stiffness  matrix  must  be  symmetric  b ^  -  b ^ ,  and  b2^ 
b^2*  The  sign  of  b2^  and  b^2  must  be  negative. 


Generally,  each  stiffness  coefficient  is  influenced  by  the  effects  of 
pile-structure  and  pile-soil  interaction.  For  example,  b^  may  be  defined 

as: 

bll  -  C1  C2 

Where  is  a  constant  depending  on  the  pinned  or  fixed  condition 
at  the  pile  head  and  C2  is  a  constant  based  on  pile-soil  interaction. 
Depending  on  the  method  used,  these  terms  may  be  calculated  separately  and 
then  multiplied  to  determine  the  pile  stiffness  coefficient,  or  the  entire 
stiffness  may  be  determined  directly. 

3.  Analysis  Method.  The  stiffness  matrix  of  each  pile  is 


transformed  from  the  local  coordinate  system  to  the  global  coordinate 
system.  All  pile  stiffness  matrices  are  then  svmsned  to  form  a  6  x  6  matrix 


representing  the  stiffness  of  the  entire  pile  group.  Applied  loads  are 
defined  as  a  set  of  three  forces  and  three  moments  acting  on  the  pile  cap. 

To  determine  displacements  of  the  pile  cap  the  following  equation  must  be 
solved: 

M  -  M  M 

Where  F  is  the  applied  load  set,  K  is  the  pile  group  stiffness  matrix,  and  U 
is  the  set  of  pile  cap  displacements.  Once  these  displacements  have  been 
determined,  the  displacements  at  the  head  of  each  pile  can  be  determined  by  a 
geometric  transformation  based  on  the  location  and  orientation  of  that  pile. 
The  following  equation  must  then  be  solved  to  determine  forces  acting  on  each 
pile  head: 

[f]  -  [b]  [u] 

Where  f  is  the  set  of  pile  loads,  b  is  the  pile  stiffness  coefficients,  and  u 
is  the  set  of  pile  head  displacements.  The  above  represents  the  basic 
analysis  of  a  pile  group.  Further  details  are  contained  in  the  user's 
manuals  for  the  various  computer  programs. 

E.  “Limitations.  Most  of  the  limitations  of  this  method  of  pile  group 
analysis  have  been  mentioned  above,  but  will  now  be  summarized.  The  method 
is  valid  for  static  analysis  of  a  linear,  elastic  system.  Applied  loads  must 
be  equivalent  static  loads,  non-linear  soil  properties  must  be  represented  by 
linear  pile  behavior.  The  other  major  limitation  is  that  the  pile  cap  is 
assumed  to  be  rigid.  Though  this  may  be  a  valid  assumption  for  a  massive 
structure,  such  as  a  dam  pier,  it  may  result  in  gross  errors  in  long,  thin 
structures,  such  as  a  U- frame  lock  monolith. 

F.  Sample  Analyses.  Several  sample  problems  are  shown  in  Appendix  E, 
solved  by  the  above  method  and  by  conventional  hand  methods. 


IV.  COMPUTER  PROGRAM  CRITERIA. 

A.  General  Requirements.  Several  related  programs  currently  use  the 
Saul  method  of  pile  group  analysis,  as  described  in  the  previous  section. 
None  of  these  programs  has  proven  completely  satisfactory  to  a  wide  range  of 
users.  The  capabilities  of  several  of  these  programs  are  shown  in  Table  1. 
The  following  paragraphs  describe  general  criteria  for  a  new  program 
utilizing  the  Saul  method.  This  program  is  intended  to  satisfy  the  widest 
possible  range  of  users.  Therefore,  it  must  incorporate  extensive 
capabilities  and  yet  maintain  a  user  oriented  format.  The  necessary 
capabilities  include  the  basic  analysis  and  comparisons  of  calculated  to 
allowable  loads.  The  required  user  oriented  features  include  convenient 
input  formats  and  user  control  over  program  operations.  A  detailed  criteria 
document  will' be  published  separately  to  fully  describe  the  required 
capabilities  of  a  new  pile  group  analysis  program. 

TABLE  1  -  Program  Capabilities  for  Pile  Group  Analysis 


New  Orleans 

LMVD 

St.  Louis 

3D 

PILE 

3D 

DOCUMENTATION 
User'  Manual 
Theorectical  Manual 
Example  Problems 

INPUT 

Data  File  Input 
Interactive  Input 
Constant  nj, 
Constant  Es 


XXX 
X  X 

XXX 

X 

XXX 
X  X 
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TABLE  1  -  Program  Capabilities  for  Pile  Group  Analysis  (continued) 


INPUT  (continued) 

Layered  £a 

Direct  bjj  Input 

Pile  Coordinate  Generation 

ANALYSIS 
Saul  2D 
Saul  3D 
Vetters 

Tension  Pile  Interation 
Checks  Calculated  vs 
Allowable  Loads 

OUTPUT 

Input  Echo 

Pile  Stiffness  Coefficients 

Pile  Group  Stiffness  Matrix 

Elastic  Center 

Structure  Deflections 

Pile  Deflections 

Pile  Forces 

Pile  Force  Components 

Sum  of  Pile  Force  Components 

Maximum  Bending  Moments  For  Pinned 

Selective  Output  Items 

GRAPHICS 

Pile  Layout 

Load  Vectors  vs  Elastic  Center 

Pile  Forces 

Pile  Load  Factors 


New  Orleans  LMVD  St.  Louis 


3D _  PILE  3D 


X  X 

X  X 

XXX 


XXX 

XXX 
XXX 
XXX 
X  X 

XXX 
X  X 

XXX 
X  X 

XXX 
Piles  X 

XXX 

XXX 
X  X 

XXX 
X  X 


B.  Program  Operation.  The  user  should  have  control  over  the  specific 
operations  to  be  performed  by  the  program  on  any  given  run.  To  provide  this 
control,  the  following  capabilities  are  required. 

1.  Timesharing.  The  program  should  run  in  the  timesharing  mode 
since  that  is  generally  more  convenient  than  batch  execution. 

2.  Input  Mode.  The  program  should  accept  interactive  input  of  data 
in  response  to  program  prompts.  It  should  also  accept  input  from  a  data 
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file.  The  interactive  input  is  useful  as  a  learning  technique  for  new  users, 
while  data  file  input  is  a  faster  method  for  experienced  users. 

3.  Output  Routing.  The  program  should  be  able  to  print  output  at 
the  timesharing  terminal  or  send  selected  data  to  a  file. 

4.  Selective  Output.  The  user  should  be  able  to  select  those  items 
he  wishes  to  see  output.  There  should  also  be  a  capability  to  print  output 
only  for  selected  piles. 

5.  Tension  Pile  Iteration.  The  program  should  be  able  to  iterate, 
at  the  user's  option,  to  account  for  the  extra  flexibility  of  piles  in 
tension. 

C.  Pile  Layout  Input.  Since  the  pile  layout  description  often 
constitutes  the  bulk  of  the  input  data,  considerable  effort  should  be  given 
to  simplifying  this  input. 

1.  Location.  Pile  locations  should  be  specified,  in  feet,  by  X,  Y, 

Z  coordinates. 

2.  Pile  Generation.  Simple  pile  generation  routines  should  be 
capable  of  easily  describing  conaon  pile  layouts  such  as  equally  spaced  piles 
between  end  points  or  rectangular  grids  of  piles. 

3.  Batter.  Batter  should  be  specified  as  a  ratio  of  vertical  to 
horizontal  distance  along  the  pile.  The  direction  of  the  batter  and  the  pile 
principal  axis  should  be  specified,  in  degrees,  as  the  angle  between  the 
batter  direction  and  the  X-axis.  Batters  and  angles  should  be  specified  in  a 
simple  manner,  such  as  specifying  a  batter  and  then  listing  all  piles  which 
have  Chat  batter. 

D.  Pile  Property  Input. 

1.  Direct  Stiffness  Input.  The  user  should  be  able  to  directly 
specify  the  coefficients  of  the  individual  pile  local  stiffness  matrix. 
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2.  Automatic  Stiffness  Calculation.  The  program  should  be  able  to 
compute  the  lateral  stiffness  coefficients  automatically  for  the  common  cases 
of  a  constant  or  linearly  varying  modulus  of  horizontal  subgrade  reaction 
(K^).  The  axial  and  torsional  stiffness  coefficients  should  be  calculated 

as  C^AE/L  and  C2JG/L,  respectively.  The  length  used  in  these 

calculations  should  be  determined  by  the  program  based  on  the  elevations  of 

the  pile  head  and  pile  tip,  and  considering  the  specified  batter  of  the  pile. 

3.  Pile  Head  Fixity.  The  lateral  stiffness  of  a  pile  depends  on  the 
degree  of  fixity  between  the  cap  and  the  pile.  The  user  should  be  able  to 
specify  this  as  fully  fixed  or  fully  pinned.  The  program  should  include  this 
fixity  when  calculating  the  lateral  stiffness  coefficients. 

E.  Pile  Allowable  Loads.  The  program  should  check  calculated  pile  loads 
against  allowables  specified  by  the  user. 

a.  Axial  Load.  The  allowable  axial  loads  specified  by  the  user 
may  depend  on  soil  capacity,  on  pile  material  capacity,  or  on  pile  buckling 
and  should  be  compared  directly  to  the  calculated  loads. 

b.  Bending  and  Axial  Load.  The  user  must  specify  allowable 
moments  about  both  principal  axes  and  an  allowable  axial  load  to  be  used  in  a 
combined  stress  equation. 

c.  Maximum  Bending  Moments.  The  program  must  calculate  maximum 
bending  moments  about  both  axes  for  use  in  the  combined  stress  equation.  The 
maximum  moments  often  occur  at  points  other  than  the  pile  head. 

d.  Overstress  Factors.  The  program  should  accept  different 
allowable  loads  for  different  load  cases  to  account  for  Group  II  loads. 

P.  Applied  Loads.  The  user  must  define  load  cases  as  sets  of  three 


forces  and  three  moments,  referenced  to  the  global  coordinate  system 


G.  Output .  The  program  should  output,  at  the  user's  option,  echoes  of 
the  pile  locations,  orientations,  properties,  and  allowable  loads;  tables  of 
calculated  pile  forces  and  combined  stress  factors  for  all  piles,  for 
selected  piles,  or  for  overstressed  piles;  and  deflections  of  the  pile  cap 
and  any  specified  points  in  space. 

H.  Pre-Processors  and  Post-Processors. 

1.  Individual  Pile  Behavior.  A  program  should  be  available  to 
calculate  pile  axial  and  lateral  stiffnesses  for  any  possible  combination  of 
pile  and  soil  properties.  This  program  should  also  be  able  to  calculate  and 
display  the  values  of  shear,  moment,  deflection  and  soil  pressure  along  the 
entire  length  of  any  pile  for  specified  pile  head  loads. 

2.  Graphical  Displays.  A  program  should  be  available  to  display  the 
specified  pile  layout,  including  batters.  It  should  alBO  be  able  to  display 
calculated  pile  forces  and  combined  stress  factors  superimposed  on  a  pile 
layout. 

3.  Pile  Interference.  A  program  should  be  available  to  check 
clearances  between  specified  piles  with  different  locations,  batters  and 
batter  directions. 

4.  Base  Slab  Analysis.  A  post-processing  program  should  be 
available  to  use  pile  forces,  transformed  to  global  coordinates,  to  help 
calculate  shears  and  moments  in  portions  of  the  pile  cap. 

V.  ADVANCED  METHODS  OF  PILE  DESIGN. 

A.  PILEOPT  Program.  PILEOPT  is  a  computer  program  intended  to  help 


determine  the  most  economical  pile  layout  possible  for  a  given  set  of  applied 
loads  and  within  constraints  specified  by  the  user.  It  was  developed  by  Dr. 
James  L.  Hill,  under  a  contract  with  the  Corps  of  Engineers.  The  program 


use a  the  same  analysis  method  described  previously  to  determine  pile  forces 
for  a  given  layout  and  applied  loads.  If  the  pile  forces  are  less  than  the 
specified  allowables,  the  program  deletes  some  piles  from  the  previous  layout 
and  reanalyzes.  The  program  also  attempts  to  choose  the  optimum  batter  for 
each  pile  group.  The  CASE  Task  Group  on  Pile  Foundations  will  furnish  a  more 
detailed  report  on  PILEOPT  at  some  future  date. 

B.  Flexible  Base  Analysis.  The  pile  analysis  method  described  above 
assumes  that  the  pile  cap,  or  structure  base  slab,  is  rigid  in  comparison  to 
the  stiffness  of  the  piles.  For  many  structures,  such  as  D- frame  lock 
monoliths,  this  is  not  a  valid  assumption,  and  the  flexibility  of  the  base 
slab  should  be  considered.  This  requires  use  of  large  programs  like  SAP  or 
STRUDL  which  can  represent  the  stiffness  of  the  structure  and  the  piles.  The 
pile  element  used  in  the  rigid  base  method  has  been  added  to  several  versions 
of  the  SAP  program  and  to  a  version  of  STRUDL .  Flexible  base  analyses  have 
already  been  performed  for  pile  founded  structures  designed  by  the  Corps  of 
Engineers.  A  more  detailed  report  on  flexible  base  analysis  will  be 
furnished  at  some  future  date. 

C.  Mon-Linear  Analysis.  One  of  the  assumptions  made  in  the  rigid  base 
analysis  method  is  that  a  pile  can  be  represented  by  a  set  of  linear 
stiffnesses.  The  actual  behavior  of  the  pile-soil  system  may  be  highly 
non-linear.  Some  existing  programs  are  capable  of  non-linear  analysis  of  a 
structure  which  is  supported  by  only  a  few  piles.  However,  for  large 
structures  supported  by  many  piles,  non-linear  analysis  is  not  currently 
practical.  A  more  detailed  report  on  non-linear  analysis  will  be  furnished 
at  some  future  date. 
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APPENDIX  A 

PILE  TYPES  AND  ALLOWABLE  STRESSES 

I.  GENERAL.  Representative  values  of  allowable  stresses  for  steel,  concrete 
and  timber  piles  are  presented  in  this  Appendix.  This  information  is 
compiled  from  data  published  by  technical  societies,  voluntary  standards 
organizations,  structural  codes,  and  Corps  of  Engineers'  guidance,  and  is 
intended  only  for  general  guidance. 

II.  TIMBER  PILES.  The  trees  most  commonly  used  for  piles  in  the  United 
States  are  Douglas  Fir,  Southern  Yellow  Pine,  Red  Pine,  and  Oak.  Timber 
piles  are  generally  the  most  economical  type  for  light  to  moderate  loads. 

They  are  available  in  lengths  from  30  to  60  ft.*  Timber  piles,  however,  are 
vulnerable  to  damage  from  hard  driving  and  to  deterioration  caused  by  decay, 
insect  attack,  marine  borer  attack,  and  abrasive  wear.  Timber  piles  are 
commonly  used  for  dolphins  and  fenders  for  the  protection  of  wharves  and 
piers  because  of  their  resilience  and  ease  of  replacement. 

A.  Allowable  Design  Stresses.  Representative  allowable  stresses  for 
pressure  treated  round  timber  piles  for  normal  load  duration  are  shown  in 
TABLE  A-l.  These  allowable  stress  values  were  derived  by  equations  specified 
by  ASTM  D2899.  "Standard  Method  for  Establishing  Design  Stresses  for  Round 
Timber  Piles".  ASTM  D2899  does  not  provide  a  method  for  establishing  the 
allowable  tensile  stress  parallel  to  the  grain.  However,  an  allowable 
tensile  stress  equal  to  the  allowable  bending  stress  may  be  used. 


*  A  table  of  factors  for  converting  inch-pound  units  of  measurement  to 
metric  (SI)  units  is  presented  on  page  3. 
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Table  A-l  Allowable  Unit  Stresses  for  Fully 
Supported,  Pressure  Treated,  Round  Timber 
Piles  -  Normal  Load  Duration  (4)  (7) 


Species 

Compression 
Parallel  to 
Grain  (psi) 
Fa(5)  (6) 

Bending 

(psi) 

Fb(6) 

Horizontal 

Shear 

(psi) 

Compression 
Perpendicular 
to  Grain 
(psi) 

Modulus 

of 

Elasticity 

(psi) 

Pacific  Coast 
Douglas  Fir 

(1)  1050 

2050 

115 

230 

1,500,000 

Southern  Pine 
(1)  (2) 

1000 

2000 

110 

250 

1,500,000 

Red  Oak  (3) 

900 

2050 

135 

350 

1,280,000 

Red  Pine 

750 

1600 

85 

155 

1,280,000 

(1)  The  working  stresses  for  compression  parallel  to  grain  in  Douglas  Fir 
and  Southern  Pine  may  be  increased  0.2  percent  for  each  foot  of  length  from 
the  top  of  the  pile  to  the  critical  section.  For  compression  parallel  to 
grain,  an  increase  of  2.5  psi  per  foot  of  length  is  recommended. 

(2)  Values  are  weighted  averages  for  longleaf,  slash,  loblolly,  and 
shortleaf . 

(3)  Values  are  weighted  averages  for  Northern  and  Southern  Red  Oak. 

(4)  The  working  stresses  in  this  table,  excepted  for  modulus  of  elasticity, 
have  been  adjusted  to  compensate  for  strength  reductions  due  to  conditioning 
prior  to  treatment.  There  piles  are  air  dried  or  kiln  dried  before  pressure 
treatment,  or  where  untreated  piles  are  to  be  used,  the  above  working 
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stresses  shell  be  increased  by  dividing  the  tabulated  values  by  the  following 
factors: 

Pacific  Coast  Douglas  Fir,  Red  Oak,  Red  Pine:  0.90 
Southern  Yellow  Pine:  0.85 

(5)  For  allowable  compressive  stresses  within  the  unsupported  length  of 
timber  piles,  see  paragraph  l.B. 

(6)  The  allowable  stresses  for  compression  parallel  to  the  grain  and 
bending,  derived  in  accordance  with  ASTM  D2899,  are  reduced  by  a  safety 
factor  of  1.2  in  order  to  comply  with  the  general  intent  of  paragraph  13.1  of 
ASTM  D2899. 

(7)  For  hydraulic  structures  the  values  in  this  Table,  except  for  modulus  of 
elasticity,  should  be  reduced  by  dividing  by  a  factor  of  1.2.  This 
additional  reduction  recognizes  the  difference  in  loading  effects  between  the 
ASTM  normal  load  duration  and  the  longer  load  duration  typical  of  hydraulic 
structures,  and  the  uncertainties  regarding  strength  reduction  due  to 
conditioning  processes  prior  to  treatment. 

B.  Allowable  Compressive  Stresses  for  Pnaupported  Piles.  The 
allowable  compressive  stress  for  cross  sections  within  the  unsupported  length 
of  timber  piles  may  be  determined  by  the  formula: 


F'a  *  rr2BS 

4.0  CKL/rA)2 


where: 


Dg  ■  pile  diameter  at  large  end  (point  of  connection  to 
superstructure)  (inches) 

“  pile  diameter  at  the  location  where  pile  is  supported  by  soil 
(inches) 

F'  “  allowable  unit  stress  in  compression  parallel  to  the  grain 
adjusted  for  KL/r  ratio,  when  F'a  <  Fa  (psi) 

E  ■  modules  of  elasticity  of  pile  species  (psi) 

L  ■  unsupported  length  of  pile  (inches) 

r ^  ”  radius  of  gyration  of  pile,  taken  at  the  location  where  the 

pile  is  supported  by  the  soil  (inches) 

K  ■  .7  for  pinned-fixed  end  conditions 

K  “  .5  for  fixed-fixed  end  conditions 

The  above  formula  is  applicable  for  a  pile  fixed  below  the  ground  surface  and 
fixed  (K  *  .5)  or  hinged  (K  “  .7)  at  the  pile  cap.  The  formula  has  a  safety 
factor  equal  to  4.0.  If  translation  of  the  pile  caps  needs  to  be  considered, 
a  critical  pile  buckling  load  may  be  determined  by  methods  outlined  in 
reference  (7)  or  by  using  the  computer  program  discussed  in  Appendix  C. 

C.  Combined  Axial  Load  and  Bending.  For  combined  axial  load  and 

bending,  stresses  should  be  so  proportioned  that: 

f./F.  ♦  V*b  *  1-0 

where: 

fa  ■  computed  axial  stress  (psi) 
fp  ■  computed  bending  stress  (psi) 

F*  “  allowable  axial  stress  (psi) 

F|>  *  allowable  bending  stress  (psi) 


The  above  formula  is  applicable  for 


a  - 


15ES 


For  > 

Da 


/n 

V  Fa 


15E8  ,  the  combined  axial  load  and  bending 


atreaa  should  be  proportional  that! 
♦  *b 

?»■.  d-fayfbi  Fb 


1.0 


where: 

F’a  is  as  defined  for  unsupported  piles 
Since  timber  pilea  are  tapered,  the  critical  section  or  point  of  maximum 
stress  may  be  at  the  tip  for  end  bearing  piles;  or  in  the  upper  region  where 
subject  to  bending,  axial  load  and  buckling;  or  at  some  point  between  for 
friction  piles. 


111.  STEEL  PILES. 

Steel  piles  in  general  are  available  in  long  lengths;  are  able  to 
withstand  hard  driving  and  penetrate  dense  strata;  and  can  carry  moderate  to 
heavy  loads.  Embedded  steel  piles  may  be  subject  to  deterioration;  by 
rusting  above  and  slightly  below  the  ground  line,  especially  in  or  near  salt 
water;  by  corrosion  if  the  surrounding  foundation  material  is  coal,  alkaline 
soils,  cinder  fills  or  wastes  from  mines  or  manufacturing  plants;  or  by  local 
electrolytic  action. 


AS 


A.  H  Piles.  H  piles  are  nondisplacement  piles  which  cause  little 


disturbance  to  the  surrounding  soil  during  driving.  H  piles  can  carry  loads 

/ 

s 

up  to  200  tons,  however,  the  usual  range  is  from  40  to  120  tons.  Their 
length,  although  basically  unlimited,  typically  ranges  between  40  to  100 
feet.  H  piles  are  easy  to  splice. 

B.  Open-End  Pipe.  Open-end  pipe  piles  can  also  be  considered 
nondisplacement  piles,  provided  they  are  augered  or  otherwise  cleaned  out  as 
they  are  driven.  They  can  be  installed  in  unliaited  lengths  and  can  carry 
moderate  loads. 

C.  Closed-End  Pipe.  Closed-end  pipe  piles  are  displacement  piles  used 
when  it  is  desirable  to  add  volume  to  and  compact  the  surrounding  soil  in 
order  to  increase  the  skin  friction  on  the  pile.  This  type  of  pile  may  cause 
heave  of  the  surrounding  piles  and  soil. 

D.  Allowable  Design  Stresses.  Allowable  design  stresses  for  steel 
piles  are  shown  in  TABLE  A-2.  Allowable  compressive  stresses  are  given  for 
both  the  lower  and  upper  regions  of  the  pile.  Since  the  lower  region  of  the 
pile  is  subject  to  damage  during  driving,  the  allowable  compressive  stress 
should  be  .28  Fy  (10,000)  psi.  This  value  may  be  increased  for  pipe  piles 
that  are  inspected  for  damage  after  driving.  Bending  and  buckling  effects 
are  usually  minimal  in  the  lower  region  of  the  pile  and  need  not  be 
considered.  The  upper  region  of  the  pile  may  be  subject  to  the  effects  of 
bending  and  buckling  as  well  as  axial  load.  Since  this  region  (from  about  IS 
feet  below  the  ground  surface  to  the  pile  cap)  is  not  usually  damaged  during 
driving,  a  higher  allowable  compressive  stress  is  permitted.  The  upper 
region  of  the  pile  is  actually  designed  in  the  same  manner  as  a  steel  column, 
with  due  consideration  to  lateral  support  conditions  and  combined  stresses. 
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TABLE  A- 2  -  Allowable  Design  Stresses  for  Steel  Piles 
Increase  Allowable  Stresses  33Z  (.67  Fy  max)  for 
Temporary  Loads.  (3) 


Code 

or 

Reference 

Compression 
at  Pile  Tip 
(psi) 

F. 

Comp.  Upper  Region 
Subject  to  Combined 
Stresses  (psi) 

Fa  (4) 

Tension 

(psi) 

Bending 

(psi) 

Fb(5)(6) 

AISC 

0.60  Fy 

0.60  Fy 

0.60  Fy 

(22,000) 

(22,000) 

(22,000) 

AASHTO 

0.472  Fy 

0.55  Fy 

0.55  Fy 

(17,000) 

(20,000) 

(20,000) 

Draft  Pile  (2) 

0.28  Fy 

0.47  Fy 

0.60  Fy 

EM 

(10,000) 

(17,000) 

(22,000) 

Recommended  (1) 

0.28  Fy 

0.47  Fy 

0.50  Fy 

0.50  Fy 

for  Hydr.  Structures 

(10,000) 

(17,000) 

(18,000) 

(18,000) 

(1)  The  recomsiended  allowable  stresses  for  hydraulic  structures  are  S/6  of 
AISC  values.  Fa  for  the  upper  region  is  based  on  an  avereage  safety 
factor  rather  than  the  variable  safety  factor  specified  by  AISC. 

(2)  Note  inconsistencies  in  Draft  Pile  EM  values  for  head  compression  and 
bending. 

(3)  Values  given  in  parenthesis  are  for  A  36  steel. 

(4)  For  allowable  compressive  stresses  within  the  unsuppported  length  of 
steel  piles,  see  paragraph  2K. 

(5)  For  combined  axial  load  and  bending,  see  paragraph  2.F. 

(6)  The  allowable  bending  stress  values  given  assume  the  compression  flange 
is  adequately  supported.  For  other  conditions  refer  to  the  allowable 
bending  stress  formulas  give  in  EM  1110-1-2101. 
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E.  Unsupported  Pile*.  The  allowable  compressive  stress 
for  an  unsupported  steel  pile,  where  Cc  >  KL  may  be  determined  by  the 


formula: 


F-  -  F« 


1  -  (EL/r)2  F» 


F.S.  i  4TT'i  E 

or  when  Cc  <  KL  by  the  formula: 


where: 


and  where: 


and  where: 


C-  -  /2t E 


allowable  axial  compressive  stress  (psi) 

specified  minima  yield  stress  (psi) 

modulus  of  elasticity  (29,000,000  psi) 

actual  unbraced  length  of  pile  from  the  pile  cap  to  the 

point  of  fixity  below  the  ground  surface  (inches) 

effective  length  factor  as  defined  by  AISC 

least  radius  of  gyration  (inches) 


F.S.  "  Factor  of  Safety 

■  Varies  from  1.67  to  1.92  for  AISC 
-  2.12  for  AASHTO 

»  2.15  for  recommended  value  for  hydraulic  structures 
F.  Combined  Axial  Load  and  Bending.  Steel  piles  subject  to  axial  load 
and  bending  shall  be  proportioned  to  satisfy  the  following  requirements: 

fa  „  °mx  gbx  „  ^y  <by  $  l  Q 
F.  U-f./F’„)Fbx  (l-f,/F'ey)Fby 


f*  +  ^bx  +  hy  £  |,o  (when  f.  $  1.5) 
*a  Fbx  rby  "?t 


mm 


where: 


and: 


it2E 


F-s-  <*b  Vrbr 


fa 

fbx  or  fby  * 


F.  “ 

Fbx  or  Fby  - 

E  - 
Lb  “ 


computed  axial  stress  (psi) 

computed  compressive  bending  stress 
about  the  x  axis  and  y  axis,  respectively 
(psi) 

allowable  axial  stress  (psi) 

allowable  compressive  bending  stressabout 
the  x  and  y  axis,  respectively  (psi) 
modulus  of  elasticity  (29,000,000  psi) 
actual  unbraced  length  of  pile  in  the  plane 
of  bending  (inches) 

effective  length  factor  as  defined  by  AISC 
in  the  plane  of  bending  (inches) 
radius  of  gyration  in  the  plane  of  bending 
(inches) 

coefficient  about  x  and  y  axes, 
respectively,  as  defined  by  AISC 
Factor  of  Safety  (see  paragraph  E) 


C.  Splices .  Splices  should  be  designed  to  develop  the  full  strength 
of  the  pile  in  compression,  tension,  and  flexure. 

IV.  CONCRETE  FILLED  STEEL  FILES. 

A.  Open-  and  Closed-End  Pipe.  Pipe  piles,  both  open-  and  closed-end, 
can  be  filled  with  concrete  to  increase  their  structural  load-carrying 
capacity.  Loads  up  to  300  tons  can  be  carried  with  this  type  of  pile. 

B.  Drilled  in  Caissons.  Drilled  in  caissons  are  typically  open-end 
pipe  piles  of  24-inch  or  30-inch  diasMter,  drilled  into  rock.  They  can  carry 
loads  up  to  300  tons.  If  an  H  Pile  core  section  is  also  used,  the  load 
carrying  capacity  can  be  increased  considerably. 

C.  Allowable  Design  Stresses.  Allowable  design  stresses  for  concrete 
filled  steel  piles  should  follow  steel  and  concrete  allowable  stresses 
specified  in  paragraphs  III  and  V  of  this  Appendix. 


D.  Pnsupported  Piles  and  Combined  Stregiet.  For  theae  conditions , 
follow  the  provisions  for  concrete  piles. 

V.  CONCRETE  PILES. 

A.  Precast  Concrete  Piles.  This  general  classification  covers  both 
conventionally  reinforced  concrete  piles  and  prestressed  concrete  piles. 

Both  types  can  be  formed  by  casting,  spinning,  or  extrusion  methods,  and  are 
made  in  various  cross  section  shapes  such  as  square,  octagonal,  and  round. 
Precast  concrete  piles  must  be  designed  and  manufactured  to  withstand 
handling  and  driving  stresses  in  addition  to  service  loads. 

1.  Conventionally  reinforced  concrete  piles  are  constructed  of 
reinforced  concrete  with  internal  reinforcement  consisting  of  a  cage  made  up 
of  several  longitudinal  bars  and  lateral  ties  of  hoops  or  spirals. 

2.  Prestressed  concrete  piles  are  constructed  using  steel  rods, 
strands,  or  wires  under  tension  to  replace  the  longitudinal  steel  used  in  the 
construction  of  conventionally  reinforced  concrete  piles.  The  prestressing 
steel  is  enclosed  in  a  conventional  spiral.  Such  piles  can  usually  be  made 
lighter  and  longer  than  normally  reinforced  concrete  piles  for  the  same 
rigidity  and  bending  strength.  Other  advantages  of  prestressed  piles  are: 

a.  Durability 

b.  Crack  free  during  handling  and  driving 

c.  High  load-carrying  capacity 

d.  High  moment  capacity 

e.  Excellent  combined  load-moment  capacity 

f.  Ability  to  take  uplift  (tension) 

g.  Ease  of  handling,  transporting,  and  driving 

h.  Economy 
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i.  Ability  to  take  hard  driving  and  to  penetrate 


hard  atrata 

j.  High  column  strength 

k.  Readily  spliced  and  connected 

B.  Cast-in-Place  Concrete  Piles.  In  general,  cast-in-  place  concrete 
piles  are  installed  by  placing  concrete  in  a  preformed  hole  in  the  ground  to 
the  required  depth.  Depending  on  foundation  conditions,  the  hole  is  usually 
lined  with  a  steel  casing  which  is  left  in  place  or  nay  be  pulled  as  concrete 
is  placed.  Since  the  concrete  is  not  subjected  to  driving  stresses,  only  the 
stress  frost  service  loads  need  be  considered  in  the  design.  Basic  types 
include  the  following:  Cased  driven  shell,  drilled-in-caisson, 

dropped- in- she 11 ,  uncased,  compacted,  auger  grouted  injected, 
caat-in-drilled-hole,  and  composite  concrete  piles.  Detailed  descriptions  of 
each  of  these  types  are  covered  in  Chapter  1  of  Reference  1. 

C.  Allowable  Design  Stresses,  the  allowable  design  stresses 
determined  in  accordance  with  the  recommended  formulas  in  this  section  relate 
to  the  structural  capacity  of  the  pile  with  an  applied  factor  of  safety.  The 
design  stresses  reflect  a  minimum  safety  factor  of  2.2  (based  on  strength 
design)  and  include  an  accidental  eccentricity  factor  of  5  percent. 

Allowable  design  stresses  for  concrete  piles  are  shown  in  TABLE  A4.  For  bond 
and  shear  allowables,  see  the  provisions  of  ACI  318-77. 


TABLE  A- 3  -  Allowable  Design  Stresses  for  Concrete  Piles 


Allowable  Stresses*  (psi) _ 

Permanent  Loads  Hydraulic  Structures 


Concrete 

Compression 


Confined** 

.40  f'c 

.35  f'c 

Unconfined 

.33  f'c 

.33  f'c 

Tension 

Plain  and  Reinforced 

0 

0 

Prestressed 

3  JFc  (250  max) 

3  JF c 

Bending  Compression 

All  Types 

.45  f'c 

.35  f'c 

Bending  Tension 

Plain 

0 

0 

Reinforced 

0 

0 

Prestressed 

3  JFZ  (250  max) 

3  JFl 

Reinforcing  Steel 

Grade  40,  50 

20,000 

20,000 

Grade  60 

24,000 

20,000 

(250  max) 


(250  max) 


♦Reduce  allowable  stresses  10Z  for  trestle  piles  and  for  piles 
supporting  piers,  docks,  and  other  marine  structures. 

♦♦Provided  the  steel  shell  confining  the  concrete  is  not  greater  than 
seventeen  inches  in  diameter;  is  fourteen  gage  (U.S.  Standard)  or 
thicker;  is  seamless  or  has  spirally  welded  seams;  has  a  yield 
strength  of  30,000  psi  or  greater;  is  not  exposed  to  a  detrimental 
corrosive  environment;  and  is  not  designed  to  carry  a  portion  of 
the  pile  working  load. 


D.  Combined  Axial  Load  and  Bending.  Por  combined  axial  load  and 

bending,  the  concrete  stresses  should  be  so  proportioned  that: 

1.  Axial  compression  and  bending: 

For  all  piles:  fa  ♦  fb  ^  1.0 
Fa  Fb 

For  prestressed  piles:  fa  +  fb  +  fpc  3  0.45  f'c 

(.35  f'c  for  hydraulic 
structures,  compression) 
fa  -  fb  ♦  fpc  *  0  (tension) 


2.  Axial  tension  and  bending: 

For  preatreased  piles:  -  fa  -  fb  ♦  fpc  fc  0 

vhere : 

fa  *  actual  axial  stress 

fb  “  actual  bending  stress 

fpc  *  effective  prestress  after  losses 

Fa  “  allowable  axial  stresss 

Fb  “  allowable  bending  stress 

3.  When  the  pile  is  designed  for  combined  axial  load  and  bending, 
the  working  stress  design  should  be  checked  using  strength  design  methods  to 
insure  that  the  required  minimum  factor  of  safety  is  achieved  in  accordance 
with  ACI  318-77. 

E.  Allowable  Design  Loads. 

1.  Laterally  Supported  Piles.  The  allowable  compressive  design 
loads  on  laterally  supported  concrete  piles  may  be  determined  by  using  Table 
A-4. 

2,  Unsupported  Piles.  Where  the  pile  extends  above  the  ground  or 
where  scour  is  expected,  the  allowable  load  must  be  reduced.  For  1/r  ratios 
up  to  120,  the  allowable  load  for  the  unsupported  pile  length  may  be 
determined  by  applying  a  reduction  factor  R,  to  the  allowable  load  for  a 
fully  supported  pile,  where  R  “  1.23  -  0.008  (1/r)  *  1.0.  If  1/r  exceeds 
120,  the  pile  should  be  investigated  for  elastic  stability.  The  effective 
pile  length  (1)  is  determined  by  multiplying  the  structural  pile  length  (L) 
by  the  appropriate  value  of  the  coefficient  K  listed  below: 

VALUES  FOR  K  FOR  VARIOUS  HEAD  AMD  END  CONDITIONS 


Head 

End  Conditions 

Condition 

Both  fixed  One  fixed 

Both  hinged 

Non¬ 

translating 

0.6  0.8 

1.0 
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F.  Other  Considerations. 

The  pile  foundation  design  should  include  other  considerations  to 
ensure  that  piles  are  installed  satisfactorily.  Some  of  these  considerations 
are  as  follows: 

1.  Pile  Dimensions.  It  is  reconmended  that  the  minimum 
dimension  be  10  inches. 

2.  Pile  Shells.  Pile  shells  or  casing  should  be  of 
adequate  stength  and  thickness  to  withstand  the  driving  stresses  and  maintain 
the  cross  section  of  the  driven  pile. 
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APPENDIX  B 

PILE  INSTALLATION  METHODS 

I.  DRIVING  BY  IMPACT  METHODS. 

Most  piles  are  installed  by  driving  with  impact  hammers.  These  hammers 
are  usually  powered  by  steam,  air,  or  diesel.  The  pile  driving  equipment 
used  should  be  adequate  to  satisfactorily  install  the  pile  to  the  penetration 
or  resistance  required  without  damage  to  the  pile.  Hammer  types  can  be 
classified  as  gravity  or  drop  hammers,  single  acting  hammers,  double  acting 
hammers,  or  diesel  hammers.  Gravity  or  drop  hammers  are  seldom  used.  They 
consist  of  a  weight  lifted  by  cable  to  a  specified  height.  The  weight  is 
released  and  the  energy,  supplied  by  gravity,  drives  the  pile.  The  single 
acting  hammer  operates  in  the  same  fashion,  only  the  weight  is  raised  by 
steam  or  air  power.  The  steam  or  air  power  permits  the  weight  to  be  raised 
and  released  much  more  rapidly  than  by  drop  hammer.  Double  acting  steam 
hammers  employ  steam  or  air  power  to  raise  the  hammer  and  to  power  the  hammer 
on  the  downward  stroke.  Diesel  pile  hammers  get  their  energy  from  the 
compression  blow  of  a  falling  weight  and  the  reaction  to  controlled 
instantaneous  burning  and  expansion  of  fuel,  which  raises  the  hammer  for  the 
next  stroke.  In  general,  the  more  driving  energy  delivered  to  the  pile, 
without  damaging  the  pile,  the  better. 

II.  PRE-EXCAVATION  METHODS. 

Pre-excavation  methods  such  as  jetting,  preboring,  augering,  or 
spudding  are  used  when  piling  must  be  driven  through  dense  or  hard  materials 
to  bearing  at  greater  depth  or  when  it  is  necessary  to  remove  an  equivalent 
amount  of  non-compresaible  soil  before  installing  displacement  piles  such  as 
closed-end  pipe,  concrete,  or  timber.  Pre-excavation  methods  will  also 


minimize  or  eliminate  the  vibration  cauaed  by  driving  which  may  damage 
adjacent  structures.  Pre-excavation  methods  should  be  used  with  care  in 
order  to  insure  the  desired  capacity  of  the  piles  being  installed  and  the 
capacity  of  the  piles  already  in  place;  and  to  insure  the  safety  of  nearby 
existing  structures. 

A.  Jetting.  Jetting  is  accomplished  by  pumping  water  through  pipes 
attached  to  the  side  of  the  pile  as  it  is  driven.  This  method  is  used  to 
install  piles  through  cohesionless  soil  materials  to  greater  depths.  The 
flow  of  water  reduces  skin  friction  along  the  sides  of  the  pile.  Air  jetting 
is  also  used.  The  pile  is  usually  jetted  to  within  a  few  feet  from  the  final 
elevation  and  then  driven.  Since  jetting  reduces  skin  friction,  it  should  be 
used  with  caution,  especially  for  tension  piles. 

B.  Predrilling  or  Augering.  Predrilling  is  used  to  produce  a  hole 
into  which  a  driven  pile  may  be  installed.  The  hole  may  be  used  to  penetrate 
difficult  materials  or  to  provide  accurate  location  and  alinement  of  the  pile. 

C.  Spudding.  Spudding  is  accomplished  by  driving  a  heavy  pipe 
section,  mandrel  or  H  pile  section  to  provide  a  hole  through  difficult  or 
hard  foundation  materials.  The  spud  is  pulled  and  the  pile  is  inserted  in 
the  hole  and  driven  to  the  required  resistance. 

III.  VIBRATORS. 

A.  Low  Frequency  Vibrators.  Low  frequency  vibrators  deliver  their 
energy  by  lifting  the  pile  and  driving  it  downward  on  each  cycle.  These 
operate  at  frequencies  of  5  to  35  cycles  per  second.  The  vibration  tends  to 
reduce  the  frictional  grip  of  the  soil  on  the  pile  and  the  pile  itself  is 
used  to  impact  the  soil  and  overcome  point  resistance.  This  method  has  found 
only  limited  use  in  driving  displacement  piles.  Use  in  the  installation  of 
nondisplacement  piles,  however,  is  increasing. 


B.  High-Frequency  Vibrators.  High  frequency  vibrators  operate  at  the 
natural  frequency  of  the  pile.  The  pile  itself  imparts  energy  to  displace 
the  soil  in  front  of  Che  pile  tip.  High  frequency  vibrators  operate  between 


40  and  140  cycles  per  second.  Displacement  piles  over  100  feet  long  have 
been  installed  using  this  method. 

IV.  CAST- IN-PLACE  PILES. 

This  method  consists  of  forming  a  hole  in  the  soil  and  filling  it  with 
concrete.  Cast-in-place  piles  may  be  cased  or  uncased.  Casings  (shells)  may 
be  driven  with  or  without  a  mandrel.  The  casings  are  driven  to  the  desired 
resistance  and  filled  with  concrete  or  the  casing  may  be  slowly  withdrawn  as 
the  concrete  is  poured  into  the  hole. 
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APPENDIX  C 


PILE  STIFFNESS  COEFFICIENTS 


I.  GENERAL. 

The  ability  o£  a  pile  foundation  to  resist  applied  loads  depends  on  the 
complex  interaction  of  the  pile  with  the  surrounding  soil.  The  nuaerous 
factors  which  affect  the  response  of  a  pile  foundation  must  be  reduced  to  a 
mathematical  representation  so  that  a  reasonably  accurate  analytical 
evaluation  can  be  performed.  The  most  common  method  of  accomplishing  this 
representation  is  to  replace  the  soil  and  pile  vith  springs  at  the 
pile-structure  interface.  Once  the  various  properties  of  the  soil-pile 
foundation  are  represented  by  equivalent  spring  constants,  it  is  relatively 
easy  to  determine  the  pile  forces  by  use  of  any  of  several  computer  programs 
currently  available.  The  difficulty  arises  in  establishing  the  equivalent 
soil-pile  springs  with  a  reasonable  degree  of  accuracy.  Two  approaches  are 
generally  used: 

-  Pile  load  test  values:  determined  by  actual  full  scale  pile  load 
tests  at  the  construction  site  or  a  nearby  site  with  similar  soil  conditions. 

-  Semi-Empirical:  determined  in  two  ways,  by  formulae  or  by  computer 
solution.  If  the  soil  modulus*  can  be  assumed  to  be  constant  or  to  vary 
linearly  with  depth,  the  equivalent  springs  can  be  determined  directly  by 
formulae  shown  in  Table  C-2.  For  more  complex  soil  systems,  a  computer 
solution  can  be  used  to  account  for  multi-layered  soils,  non-linear  variation 
of  soil  modulus,  and  inelastic  soil  behavior  by  analyzing  a  single,  isolated 
pile  using  known  soil  parameters  for  the  site. 

♦NOTE:  In  this  Appendix  Es  refers  to  the  horizontal  soil  modulus. 


The  equivalent  springs  or  stiffness  coefficients,  determined  by  the 
methods  described  above  are  the  “b"  terms  of  the  pile  stiffness  matrix  as 
described  in  paragraph  III  D2  of  the  text.  Generally,  these  terms  can  be 


defined  as: 


bij  "  C1C2 


where  Cj  is  a  constant  which  depends  on  the  fixity  of  the  pile  head  to  the 
structure.  For  most  applications,  a  fixity  condition  of  fully  pinned  or 
fully  fixed  is  assused.  is  a  constant  based  on  the  pile-soil 
interaction  and  is  determined  by  one  of  the  methods  mentioned  above. 

Values  for  the  fixity  constant  Cp  for  soils  with  a  constant  or  a 
linearly  varying  soil  modulus,  Eg,  are  shown  in  Table  C-l.  A  theoretical 
derivation  of  these  values  can  be  found  in  references  16  and  17. 


TABLE  C-l 

Pile  Fixity  Constants,  For  Soils  With 
Constant  or  Linear  Variation  of  Soil  Modulus 


Pile  Stiff.  Coeff. 


Values  of  Ci 


Const.  Es 
Pile  Head  Fixitj 


Linear  Var.  Es 
Pile  Head  Fixit? 


1.075  0.411 


1.075  0.411 


The  pile  stiffness  coefficients  and  representing  the 
axial  and  torsional  stiffness,  respectively,  of  the  pile  are  not  shown  in 
Table  C-l.  These  two  coefficients  are  assumed  to  be  not  affected  by  the  pile 
head  fixity  and,  therefore,  are  not  shown.  For  additional  discussion  of  the 
axial  and  torsional  pile  stiffness  coefficients,  see  sections  III  and  IV, 
respectively,  of  this  appendix. 

The  pile  stiffness  coefficients  can  be  affected  by  many  factors  other 
than  the  pile  fixity  constants  (Cj)  and  the  pile-soil  stiffness  constants 
(Cj)*  The  major  factors  are  mentioned  here  but  a  detailed  discussion  is 
beyond  the  scope  of  this  appendix.  The  following  items  could  influence  the 
pile  stiffness  coefficients: 

-  Group  effect:  close  spacing  of  piles  in  a  large  group  can  reduce 
the  lateral  capacity  for  the  group. 

-  Position  in  group:  a  pile  may  exhibit  different  stiffness 
depending  on  its  location  in  the  group. 

-  Cyclic  loading:  repeated  application  of  static  loads  on  a  pile  can 
cause  greater  deflections  of  the  pile  than  the  application  of  a  sustained 
static  load  of  equal  magnitude. 

-  Vibratory  or  dynamic  loading:  statically  loaded  piles  subjected  to 
vibrations  or  dynamic  loads  may  deflect  significantly  more  than  with  the 
static  load  only. 

-  Driving  a  pile  into  closely  spaced  group:  when  piles  are  driven  in 
an  area  that  already  contains  closely  spaced  pile,  the  soil  density  within 
the  pile  group  can  be  affected. 

-  Sheet  pile  cutoff:  sheet  pile  used  to  inclose  pile  groups  may 
change  the  distribution  of  stress  in  the  soil. 
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-  Water  table  and  seepage:  the  groundwater  table  and  seepage  can 
influence  the  lateral  soil  modulus. 

-  Pile  length:  short  rigid  piles  act  differently  than  long  flexible 
piles.  This  report  assumes  piles  are  long  enough  to  act  in  a  flexural  mode 
(non-dimensional  length  L/T  is  greater  than  5,  as  defined  by  Reese  (17)). 

-  Stiffness  of  pile  cap:  the  flexibility  of  the  pile  cap  will 
influence  the  distribution  of  load  to  the  piles. 

For  additional  discussion  of  the  factors  mentioned  above,  see  reference 

18. 

The  remainder  of  this  appendix  will  deal  with  determination  of  the  pile 
stiffness  constants  (C^)  without  regard  to  the  items  briefly  referred  to 
above . 

II.  LATERAL  STIFFNESS. 

A.  General .  For  structures  which  experience  lateral  loads  of  any 
significance,  the  correct  representation  of  the  lateral  stiffness  of  the 
foundation  in  the  anlysis  is  critical.  This  representation  must  include  the 
resistance  of  the  pile  to  lateral  translation  and  rotation  and  the  coupling 
effects.  These  stiffnesses  are  inserted  in  the  pile  stiffness  matrix  as  the 
terms  bu,  b2J,  b44,  b55,  b1$,  b24,  b42,  and  b$1.  These 

terms  can  be  determined  either  by  pile  load  tests  or  by  semi-empirical 
methods . 

B.  File  Load  Tests.  The  pile  stiffness  coefficients  can  be  determined 
by  full  scale  pile  load  tests  at  the  construction  site  or  a  nearby  site  with 
similar  soil  conditions.  However,  pile  load  tests  may  not  be  practical  for 
design  for  several  reasons: 


1.  The  tests  ere  usually  very  costly  and  time  consuming  and  may 
not  be  economically  feasible  for  small  to  medium  size  jobs. 

2.  Normally,  pile  load  tests  at  the  construction  site  are  not 
conducted  until  construction  is  well  underway.  Since  pile  analysis  and 
design  must  be  accomplished  well  in  advance  of  construction,  data  obtained 
from  load  tests  could  not  be  used  for  design  but  only  for  verification  or 
modification  of  the  pile  design. 

3.  With  restricted  site  areas,  the  pile  load  tests  can  be  in  the 
way  of  other  construction  and,  in  some  instances,  actually  delay  construction. 

C.  Semi-Empirical  Methods.  These  methods  can  be  categorized  as 
analytical  (using  formulae)  or  as  numerical  (using  a  computer  solution). 

1.  Analytical  Method.  If  a  soil  system  can  reasonably  be  assumed 
to  have  a  soil  modulus  that  varies  linearly  with  depth  or  that  is  constant, 
then  the  lateral  stiffness  constants  can  be  calculated  using  prescribed 
values  or  ranges  of  values  of  the  soil  modulus.  Shown  in  Table  C-2  are 
formulae  for  calculating  the  lateral  stiffness  terms  (b^  and  b22>,  the 
rotational  stiffness  terms  (b^  and  b^) ,  and  the  coupling  stiffness 
terms  (bjj,  **24’  **42’ 

These  terms  are  defined  as: 

b^  is  the  force  required  to  displace  the  pile  head  a  unit 

distance  along  the  local  1  axis 

b22  is  the  force  required  to  displace  the  pile  head  a  unit 

distance  along  the  local  2  axis 

b-jj  is  the  force  required  to  displace  the  pile  head  a  unit 
distance  along  the  local  3  axis 
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is  the  moment  required  to  displace  the  pile  head  a  unit 
rotation  along  the  local  1  axis 

is  the  moment  required  to  displace  the  pile  head  a  unit 
rotation  along  the  local  2  axis 

*b^  is  the  force  along  the  local  1  axis  caused  by  a  unit  rotation 
of  the  pile  head  around  the  local  2  axis 
*b24  is  the  force  along  the  local  2  axis  caused  by  a  unit  rotation 
of  the  pile  head  around  the  local  1  axis 
*b^|  is  the  moment  around  the  local  2  axis  caused  by  a  unit 
displacement  of  the  pile  head  along  the  local  1  axis 
*b^2  is  the  moment  around  the  local  1  axis  caused  by  a  unit 
displacement  of  the  pile  head  along  the  local  2  axis 

♦Since  the  stiffness  matrix  must  be  symmetric  bjj  *  b^^  and 
b24  *  b42'  The  3ign  of  b2if  and  b^2  must  be  negative. 

TABLE  C-2 

PILE  STIFFNESS  COEFFICIENTS 


Pile  Stiff.  Coeff. 

Constant  Es 

Lin.  Var.  ] 

bll 

C t Fit,. 

V 

b22 

c ,  -£&  - 

Ct  Els. 

V 

b44 

C  l  . . 

C;  Ell- 
Ti 

b55 

c}  JLs— 

b15 

Cl  _£s_. 

1  2 
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TABLE  C-2  (Continued) 
PILE  STIFFNESS  COEFFICIENTS 


Pile  Stiff.  Coeff. 

Constant  Es 

Lin. 

Var.  Es 

b24 

-Eg-. 

N 

u 

1 

E  Ti 

V 

bn  2 

-  Ci.-i.a- 

1 

o 

Ell, 

V 

b51 

C, 

Els 

V 

where : 
Cl 

Tl  = 


is  the  pile  fixity  constant  as  shown  in  Table  C-l 
from  one  "b"  term  to  another. 

(in.)  ;  T2  =  5  /El 2 

v-„ 


and  varies 

(in. ) 


"h 


Es 


Pi 


E 

I 


is  the  constant  of  horizontal  subgrade  reaction  or  the  change 
in  the  soil  modulus  with  depth  (lb/in3). 

is  the  horizontal  soil  modulus  (lb/in2). 


is  the  modulus  of  elastically  of  pile  (lb/in2), 
is  the  moment  of  inertia  of  pile  (in2*). 


Subscripts  1  and  2  for  I,  T,  and  p  refer  to  the  local  pile  axes.  See 
Figure  C-l. 
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LOCAL  PILE  AXIS 


IacoI  avis  I 


local  axis  2 


FIGURE  C-I 


The  constant  of  horizontal  subgrade  reaction  (n^)  or  the  horizontal 
soil  modulus  (Ea)  can  be  obtained  using  methods  shown  in  Appendix  D.  These 
methods  are  baaed  on  work  by  Terzaghi  (19),  Broms,  and  others  and  include 
corrections  for  pile  group  effect  and  for  cyclic  loading.  These  methods  of 
soil  modulus  are  satisfactory  if  the  variation  of  the  soil  modulus  with  depth 
can  be  reasonably  approximated  as  constant  or  linear.  Many  foundation  strata 
fall  in  this  category  (12,  17)  and  can  be  conservatively  represented  by  using 
a  "bracket"  approach  to  the  pile  design.  This  means  the  pile  foundation  is 
analysed  with  weak  pile  stiffness  coefficients  and  strong  pile  stiffness 


coefficients,  where  "weak*1  end  "strong"  refer  to  the  range  of  soil  modulus 
that  could  reasonably  be  expected  for  a  particular  soil.  In  cases  where  the 
simplified  assuaptions  are  not  valid,  computer  solutions,  are  needed. 


2.  Numerical  Solution  by  Computer.  Most  analytical  methods  are 
based  on  a  pile-soil  system  similar  to  a  beam  on  elastic  foundation.  These 
methods  assume  that  the  soil  can  be  represented  by  a  series  of  closely 
spaced,  independent  springs.  The  pile-soil  relationship  can  be  expressed  by 
a  4th  order  differential  equation  which  can  be  solved  for  specific  cases  by 
making  certain  assumptions.  There  are  several  computer  programs  available 
which  can  be  used  to  determine  the  pile  stiffness  coefficients  for  a  single 
pile.  Some  of  the  most  useful  programs  are  discussed  in  the  following 
paragraphs . 

a.  "Pile  Head  Stiffness  Matrices."  This  program  was  written  by 
Dr.  William  Dawkins  for  WES.  This  program  is  intended  to  be  used  to  analyze 
a  single  pile  to  determine  the  stiffness  coefficients  for  input  to  a  general 
pile  foundation  analysis  program.  The  procedure  used  is  a  one-dimensional 
analysis  of  a  beam  on  an  elastic  foundation  where  the  soil  is  represented  as 
discrete  springs.  The  soil  springs  are  calculated  by  the  program  based  on  a 
variation  of  the  lateral  soil  modulus  of  ■  a  ♦  bz11. 

where:  Es  ■  lateral  soil  modulus 
a,  b  ■  constants 
z  *  depth  below  ground  surface 

The  values  of  "a",  "b",  "z",  and  "n"  are  input  by  the  user.  Any  degree  of 
fixity  for  the  pile  head  to  the  pile  cap  can  be  considered  with  this 
program.  Output  consists  of  the  actual  pile  stiffness  coefficients  ("b" 
terms)  and  may  be  used  directly  as  input  to  a  general  pile  foundation 
analysis  program.  Disadvantages  of  this  program  are  that  the  user  must  know 


the  variation  of  soil  modulus  with  depth  and  that  the  soil  springs  used  are 
linearly  elastic.  Also  the  current  version  of  the  program  does  not  contain 
provisions  for  variation  of  the  pile  stiffness  with  depth  or  for  the 
application  of  axial  loads.  For  information  on  this  program  see  reference  4. 

b.  "Analysis  of  Laterally  Loaded  Piles  by  Computer”.  Several 
programs  are  available  which  can  be  used  to  determine  the  pile  stiffness 
coefficients  if  the  values  of  the  soil  springs  can  be  determined  by  other 
means.  The  values  for  the  soil  springs  are  input  to  the  program  and  the 
springs  are  treated  as  completely  elastic  or  elastic-plastic,  depending  on 
the  program's  capabilities.  Any  variation  of  the  soil  modulus  with  depth  can 
be  represented  by  inputting  the  proper  values  for  the  discrete  soil  springs. 
Axial  loads  and  variation  of  pile  stiffness  with  depth  can  generally  be 
accounted  for  in  these  programs.  Output  usually  consists  of  values  for  the 
deflection,  moment,  and  soil  reaction  for  specified  increments  along  the  pile 
model.  The  pile  stiffness  coefficients  ("b"  terms)  can  be  obtained  by 
applying  displacements  and  rotations  to  the  pile  model  and  then  using  the 
output  of  forces,  moments,  and  displacements  to  determine  the  appropriate  "b" 
terms.  For  example,  to  determine  "b^",  which  is  the  force  required  to 
displace  the  pile  head  a  unit  distance  along  axis  1,  apply  a  force  at  the 
pile  head  along  axis  1.  Then  is  equal  to  the  displacement  of  the 

pile  head  divided  by  the  applied  force.  The  other  "b"  terms  can  be 
calculated  in  similar  manner. 

It  should  be  noted  that  when  the  soil  response  to  applied  loads  is 
non-linear  (as  it  is  assumed  to  be  in  this  discussion)  the  pile  head  moments 
and  displacements  will  vary  non-linearly  with  the  forces  applied  to  the 
pile.  For  example,  if  the  applied  lateral  force  along  the  pile  axis  1  is 
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increased  linearly,  Che  pile  displacements  will  increase  non-linearly  and 
therefore''bn"vill  not  be  a  constant  but  will  vary.  In  order  to  account  for 
this  non-linearity,  the  designer  should  determine  the  sensitivity  of  the 
particular  foundation  to  variations  in  applied  pile  loads.  This  can  be 
accomplished  by  comparing  results  from  the  application  of  small  and  large 
loads  in  the  single  pile  analysis.  If  the  foundation  is  determined  to  be 
sensitive  to  the  load  variations  then  the  designer  could  account  for  this  in 
the  analysis  by  using  a  bracket  approach  for  the  "b"  terms  or  by  determining 
one  set  of  "b"  terms  which  reflect  expected  applied  pile  loads. 

One  of  the  more  useful  single  pile  analysis  computer  programs  is 
"Analysis  of  Laterally  Loaded  Piles  by  Computer,"  written  by  Dr.  Lymon 
Reese.  For  this  program,  soil  properties  are  defined  by  a  set  of  curves 
which  give  soil  reaction  as  a  function  of  pile  deflection.  The  lateral 
resistance  of  the  soil  is  represented  by  non-linear,  discrete  springs  called 
p-y  curves.  These  curves  have  been  constructed  for  various  soil  conditions 
based  on  pile  tests,  theories  for  the  behavior  of  soil  under  stress,  and 
failure  mechanisms  for  pile-soil  systems.  The  program  performs  an  iterative 
solution  which  consists  of  finding  a  set  of  elastic  deflections  of  the  pile 
which  simultaneously  satisfy  the  specified  non-linear,  resistance-deformation 
relations  (p-y  curves)  of  the  soil  and  the  elastic  bending  properties  of  the 
pile.  This  program  can  account  for  changes  in  pile  types  with  depth;  applied 
axial  loads;  a  layered,  non-linear,  soil  system;  and  any  degree  of  fixity  of 
the  pile  head  to  the  structure.  The  p-y  eurves  are  a  necessary  input  to  this 
program. 
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Another  program,  "Generation  of  Soil  Resistance  versus  Pile 
Movement  Curves,"  is  available  which  will  generate  the  required  family  of  p-y 


curves  for  a  particular  type  of  soil.  Input  to  this  program  consists  of 
several  parameters  such  as  pile  diameter,  soil  stress-strain  relationship, 
unit  weight  of  soil,  internal  friction  angle  if  sand  or  cohesion  if  clay,  and 
relative  density  if  sand  or  consistency  if  clay.  Once  the  p-y  curves  are 
determined  and  input  into  the  pile  analysis  program,  the  pile  stiffness 
coefficients  can  be  determined  from  the  output.  The  output  is  in  the  form  of 
pile  deflection  and  moment,  soil  modulus,  and  soil  reaction  values  as  a 
function  of  depth  along  the  pile.  The  only  major  disadvantage  of  this 
program  is  that  the  user  must  calculate  the  pile  stiffness  coefficients  from 
the  output.  The  lateral  stiffness  coefficients,  b^  and  b^,  m*y 
calculated  using  applied  loads  and  pile  head  displacements  (computer  output); 
and  the  coupling  coefficients,  b^  and  bjp  ®«y  likewise  be  calculated 
using  the  displacements  and  moments.  Since  the  stiffness  matrix  must  be 
8yimaetric,  b^2  ■  b^  and  bjj  “  b^^.  The  calculation  of  the 
rotational  stiffness  coefficients,  b/</(  and  b^g,  is  accomplished  by 
specifying  a  slope  (rotation)  at  the  top  of  the  pile  and  then  using  this 
slope  and  the  outputted  moments  to  determine  b^  and  b^j.  Dr.  Reese's 
programs  have  been  adapted  for  the  Corps'  timesharing  library.  For 
documentation  of  these  programs  see  reference  23.  For  an  example  using  these 
programs  to  determine  the  pile  stiffness  coefficients,  see  the  sample  problem 
at  the  end  of  this  appendix. 

E.  Summary  and  Recommendations.  Calculation  of  the  lateral  pile 
stiffness  coefficients  for  use  in  a  pile  foundation  analysis  can  be 
accomplished  by  using  data  from  pile  load  tests  or  by  mathematically 
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analyzing  a  single  pile.  If  soil  parameters  are  not  well  defined,  a 
"bracket"  approach  should  be  used  for  the  analytical  method  to  account  for 
the  maserous  unknowns  and  assumptions  involved.  Results  from  the  analysis  of 
a  single  pile  using  one  of  the  computer  programs  discussed  can  be  used  to 
verify  the  validity  of  the  upper  and  lower  limits  of  the  "bracket".  Another 
approach  when  soil  parameters  are  sufficiently  defined  through  testing,  is  to 
develop  a  set  of  stiffness  coefficients  for  the  anticipated  loads  using  one 
of  the  programs  discussed  and  then  use  these  coefficients  for  the  pile 
foundation  analysis.  In  any  case,  where  large  numbers  of  pile  are  used  or 
the  subsurface  conditions  are  out  of  the  ordinary,  analytical  assumptions  and 
results  should  be  verified  and/or  modified  by  actual  pile  load  tests. 

III.  AXIAL  STIFFNESS. 

A.  General .  The  axial  stiffness  of  a  pile  depends  on  many  factors 
such  as  the  modulus  of  elasticity  of  the  pile,  the  pile  area,  the  pile 
length,  the  pile  tip  deflection,  the  distribution  of  axial  skin  friction 
along  the  pile,  and  the  percentage  of  axial  load  transmitted  to  the  tip. 

Many  of  these  factors  are  greatly  affected  by  other  related  items  such  as 
type  of  pile  hammer,  level  of  the  water  table,  soil  density,  etc.  Some  of 
the  factors  mentioned  above,  such  as  the  pile  length,  area,  and  modulus  of 
elasticity  are  easily  determined  while  some  of  the  others  are  more  difficult 
to  ascertain. 

B.  Tip  Deflection  and  Distribution  of  Axial  Forces  Along  the  Pile. 

The  pile  tip  deflection  under  load  and  the  manner  in  which  the  axial  force  in 
the  pile  is  transmitted  to  soil  are  interrelated  and  can  have  a  great  effect 
on  the  axial  stiffness  of  the  pile.  Research  has  indicated  (21  and  22)  that 
the  amount  of  load  resisted  by  skin  friction  along  the  pile  is  dependent  on 
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the  amount  of  pile  tip  deflection.  Predicting  pile  tip  deflection  accurately 
ia  very  uncertain.  Moat  of  the  group  pile  analyaea  to  date  have  aaaumed 
that,  for  compreaaive  loada,  pile  tip  deflection  under  8ervice  loada  ia 
negligible.  For  thia  assumption,  the  axial  stiffness  can  be  assumed  to  be 
AE/L  for  an  end  bearing  pile  with  no  load  resisted  by  skin  friction  and  2AE/L 
for  a  friction  pile  with  no  end  bearing  load  transfer.  These  axial 
stiffnesses  are  analogous  to  column  effective  lengths  of  L  and  L/2.  Based  on 
a  particular  load  distribution  between  end  bearing  and  skin  friction,  the 
modifier  for  the  axial  stiffness  could  vary  between  one  and  two.  It  must  be 
emphasized  that  the  above  discussion  applies  to  pile  with  compressive  loads 
and  negligible  tip  deflection.  Tests  have  shown  that  pile  having  tensile 
loads  are  less  stiff  (as  much  as  a  50  percent  reduction)  than  piles  with  a 
compressive  load.  Furthermore,  deflection  of  the  pile  tip  of  a  relatively 
small  amount  can  cause  the  axial  stiffness  to  be  significantly  different. 
Additional  research  needs  to  be  done  to  more  definitely  predict  the  axial 
pile  stiffness  for  piles  with  tip  deflection.  In  the  absence  of  better  data, 
the  values  for  axial  pile  stiffness  shown  in  Table  C-3  have  been  used  by  some 
designers. 

A  computer  program  developed  by  Drs.  Lymon  Reese  and  H.M.  Coyle  can 
be  used  to  compute  load-displacement  relationships  for  axially  loaded  piles. 
The  load  transfer  curves  used  in  the  program  to  relate  skin  friction  to  the 
axial  displacement  of  the  pile  are  based  on  semi-empirical  criteria.  For 
more  information  on  this  program  see  reference  23. 


TABLE  C-3 


AXIAL  PILE  STIFFNESS  COEFFICIENTS 
(Assuming  No  Tip  Deflection) 


Condition 

b« 

Compressive  load. 

AE 

end  bearing  pile 

L 

Compressive  load, 

2AE 

friction  pile 

L 

IV.  TORSIONAL  STIFFNESS. 

For  e  three  dimensional  pile  group,  a  torsional  pile  constant  (b^g) 
can  be  defined  which  relates  the  rotation  of  the  pile  in  a  plane 
perpendicular  to  its  longitudinal  axis  to  an  applied  torque.  This  can  be 
expressed  as 

b66  *  ct  JG 

L 


where: 

Cj  is  a  constant  which  describes  the  distribution  of  torsional 

shear  to  the  soil  and  the  transfer  of  torsional  shear 
resistance  from  the  pile  to  the  structure. 

J  is  the  polar  moment  of  inertia  of  the  pile 

6  is  the  shearing  modulus  of  elasticity 

L  is  the  length  of  pile 

Unless  the  pile  group  is  small  (say  less  than  10  piles),  the  torsional 
stiffness  of  the  individual  pile  appears  to  have  little  effect  on  the 
stiffness  of  the  pile  group  and  can  be  conservatively  assuned  to  be  xero. 


EXAMPLE  PROBLEM: 

ANALYSIS  OF  SINGLE  PILE  TO  DETERMINE  PILE  STIFFNESS  COEFFICIENTS 
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STEP  #1.  Dtteniu  the  p-y  curves  for  the  soil  streta  shown  in  Figure  C-2. 
Use  computer  program  CORPS/10004  (MAKS)  to  generate  the  p-y  curve  data.  This 
program  also  has  graphics  capabilities  which  can  be  used  to  plot  the  input 
data  and  the  p-y  curves.  Four  p-y  curves  are  generated t  one  at  the  ground 
surface,  one  at  the  bottom  of  the  sand  layer,  one  near  the  top  of  the  clay 
layer,  and  one  at  the  bottom  of  the  clay  layer.  The  plot  of  the  input  data 
for  the  example  is  shown  in  Figure  C-3,  Figures  C-4,  C-5,  and  C-6  are  plots 
with  different  scales  of  the  p-y  curves  developed  by  the  program.  Figure  C-7 
is  a  listing  of  the  output  file  for  the  example.  For  a  complete  program 
description  and  variable  definitions  see  reference  23. 


FIGURE  C-4 :  PLOT  OF  OUTPUT 
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FIGURE  C-5:  PLOT  OF  OUTPUT 
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C-6:  PLOT  OF  OUTPUT 
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FIGURE  C-7 :  LISTING  OF  OUTPUT 
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STEP  #2.  Once  the  p-y  curves  are  determined,  the  force-displacement  data  for 
the  curves  are  used  to  represent  the  soil  behavior  and  are  input  to  a  program 
to  analyze  a  single,  laterally-loaded  pile.  For  this  analysis  use  program 
CORPS/10001  (COM62).  An  additional  p-y  curve  (#0)  has  been  added  at  the  top 
of  the  pile.  This  set  of  data  merely  represents  the  absence  of  soil  at  the 
top  of  the  pile.  Program  COM62  uses  a  linear  interpolation  between 
succeeding  p-y  curves  that  are  input.  Using  the  p-y  curves  as  specified 
force-deformation  relations  for  the  soil  and  the  pile  flexural  stiffness,  the 
program  solves  the  differential  equations  representing  the  pile-soil 
relationship.  This  is  performed  in  an  iterative  fashion  until  the  elastic 
behavior  of  the  pile  matches  the  specified  soil  behavior.  Four  runs  of 
computer  program  COM62  are  required.  Two  runs  are  required  to  determine  the 
respective  pile  stiffness  coefficients  for  each  principal  axis  of  the  pile. 

In  the  first  run  only  a  lateral  force  is  applied  along  pile  axis  1  to  the  top 
of  the  pile  In  the  second  run,  a  slope  along  the  pile  axis  1  is  applied  at 
the  top  of  the  pile.  The  graphics  associated  with  program  COM62  give  plots 
of  the  input  as  well  as  the  soil  pressure,  moments,  and  deflections  for  the 
pile.  Plots  of  the  input  data  for  the  first  two  runs  are  shown  in  Figures 
C-8  and  C-9  and  a  plot  of  the  output  for  the  runs  1  and  2  are  shown  in 
Figures  C-10  and  C-ll  (respectively).  A  listing  of  selected  output  is  shown 
in  Figure  C-12.  For  a  complete  program  description  and  variable  definitions 
see  reference  23. 
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UT  DATA 


FIGURE  C-9:  PLOT  OF  INPUT  DATA 


iMauanou 


>10:  PLOT  OF  OUTPUT  FOR  RUN  1 


FIGURE  C-ll:  PLOT  OF  OUTPUT  FOR  RUN  2 


xnaino  aaioa' 


2  *4  «4  v*  *4 

S  UthlWU 

JO900 

«*4 

hNNIUIU 
111  •  •  •  • 

9000 


t  0  0  •  • 

•  999 


s 

n 

-i 

9 

i.i  «  «  »  • 

9999 

OMOUM0 


I  UlUIUIIil 

ggggg 

•BBWjJ 
cap-inn 
*  •  •  • 

2212 


§*S8gg 

K^aeoSni 

axunicucu 


•  uiuiui 

St  cgi 

f£x*  2n® 

3  ♦  ♦  •  ♦ 

O  •••<» 


RUN  2  -  SPECIFIED  SLOPE,  STRONG  AXIS 


01- 


S  UUWWWWWWWWWWWWWUU 
«099S«0909fi09S9900 
J99999999999999999 

•*0101010101010101010101010101010101 

m«  ........  ........ 

99S99999999999999 

nnnnnrmnnnnnn 

9999999999999 


\ 

m 


UllilUlilUUUI 

▼  oooioinimoo 
nitottvv*) 
r»  ▼  cu  «*  •*  •*  oi 
•«  m  u>  coco  co  in 


UiUJ 
COOl 
01 IA 
•*m 


* 

. . . 

9999999999 


UU 

or  ® 

▼  CO 
•*01 
win 


.  .  .  . 


UUI 

▼  •4 

OOl 

9UI 

in  ▼ 

.  • 


01 


s 

00 

mi 


ninnnn' 

99990991 


till 


I  I 

▼  T 
99 


cn 

U1  ••••♦-- 


UUUUUUUUUUUUU 
r-r-cooimo-oocooi^m^ni 
in9-*roc^,<ro>'r»*nt^s3C7) 
o-  9  *«  r-  9  *4  w  •*  w  w  o-  o> 
ro**oiDcs«*oicsnoioioini 


9999999 


IVI- 

A99 

mi 

I  uu 
zww 

MOO  00 
*10  IA 
cnn 
.  * 
09 
I  I 
01*4 
99 


c^o-o-o-t-n-N^p- 

999999999 

UUUUUUUUU 

in  w  w  w  r-  •*  w  o»  ▼ 

9  00  09  ^90)01^9 

co  w  co  co  ▼  oo  n  ▼  oi 

9  9  9  9  01  0090- CO 


P-P-P-P-COW 

999999 


uuu 

888 
oo  in  ▼ 
▼  moi 


999999 _ 

I  I  I  I  I  I  I  I  I 

999999999 


99  I 

I  I  I 

h99 

999 

I 

UUU 
o-  com 
•4  9  w 
co  mco 


uuu 

wr-co 

wwco 

▼  com 

•4  0-01 

... 

99 
I  I  ( 
999 
999 

UUU 

▼  9H 

▼  W*4 
9  COW 

Ol  •*  w* 


y  .  t.i  i.i  i.i  i.i  iii  i.i  i.i  i.ii.i  i»i  i.i 
OZ^OJOOWWIAWmWmOO 
MHHOOCO<4VWW'Hlt-9 
i>  ^(uiAoxotTisono 
«>HOiMntnn)HHi«w 

c  ........... 

m  99999999999999999 

O  I  I  I  I  I  I  I  I  I  I  I 

U.  m  M  *4  •*  m4  ni  01 01 01 01 01 01 01 M  01  (U 

Z  9999999999999999 

M 

uuuuuuuuuuuuuuuu 
on®on®om«Noooo®o 
com9com9com9com9com9co 
com9com9*<mwcooo9*4mwco 
•4mwcow*4*4M«4M*«oioioioicu 


M- 

DX 
O 


99999999) 


i 


C31 


STEP  3.  Using  the  output  from  first  run,  the  stiffness  coefficients  b^, 

,  and  b^  can  be  determined: 

b  ■  applied  lateral  force  along  pile  axis  1  ■  10.000  lbs.  ■  34,258  lb/in 
1  output  displacement  at  top  of  pile  0.2919  in. 

b51  *  output  moment  at  top  of  pile  *  938,500  in- lb  “  3.215xl06  in- lb/in 
output  displacement  at  top  of  pile  0.2919  in 


b 


15 


b 


51 


Using  the  output  information  from  second  run,  the  stiffness  coefficient  b^ 
can  be  determined: 


b,,  *  output  moment  at  top  of  pile  _  9.685x10**  in-lb  _  9.685x10^  in-lb 

applied  rotation  (slope)  ”  0.10  rad  rad 

The  stiffness  constants  b 22 »  b^,  b^»  and  b^  can  be  similarly 

obtained  by  making  a  third  and  fourth  run  of  COM62  with  the  lateral  force  and 

slope  applied  with  respect  to  the  other  pile  axis.  These  values  for  the 

sample  problem  are: 

b22  -  16,113  lb/in 

b„ ,  -  1.385xl06  in-lb 

24 

b,.  -  4 . 236xl07  in-lb 

44  - 2 

rad 

b42  -  b42 

The  stiffness  constants  b^  (axial)  and  b^^  (torsional)  are  calculated  by 
other  means  (see  sections  III  and  IV  of  this  appendix) . 
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APPENDIX  D 


SOIL  MODULUS  FOR  LATERALLY  LOADED  PILES 

The  leteral  pile  stiffness  coefficients  discussed  in  Appendix  C  can  be 
directly  coaputed  when  E#  is  a  constant  or  linearly  varying  in  the  aedia  in 
which  the  pile  is  eabedded.  It  is  generally  assisted  that  for  homogeneous 
cohesive  soils  EB  is  constant  and  for  homogeneous  cohesionless  soils  E# 
varies  linearly  with  depth.  In  this  Appendix,  typical  values  for  computing 
E(  for  homogeneous  cohesive  and  cohesionless  soils  are  provided.  The 
structural  engineer  must  rely  on  the  geotechnical  engineer  to  obtain  the 
values  of  the  soil  modulus  and  the  coefficient  of  horizontal  subgrade 
reaction. 


Definitions  and  Nomenclature 


*1  -  Coefficient  of  horizontal  subgrade  reaction  (lbs/ft^)  (ratio  of 

pressure  (lbs/ft2)  at  a  point  to  the  displacement  (ft)  at  the 
point)  for  a  1  foot  wide  pile  embedded  in  clay. 

K||  -  Coefficient  of  horisontal  subgrade  reaction  (lbs/ft^)  (Kh  * 

Ea  -  Soil  modulus  (lbs/ ft2)  -  ratio  of  soil  resistance  (p)  (lbs/ft) 

to  pile  movement  (y)  (ft),  Ea  “  K^B). 

n||  -  Constant  of  horisontal  subgrade  reaction  (lbs/ft^)  for  a  pile  1 

foot  wide  embedded  in  sand  (K|,  *  n^Z/B,  Ea  ■  n^s). 

B  -  Width  of  pile  (ft). 

Z  -  Depth  below  ground  surface  (ft). 

qtt  -  Unconfined  compressive  strength  of  clay  (lbs/ft2). 

*1  *  *2  -  Reduction  factors 

f  -  Unit  weight  of  soil 


Homogeneous  Soil* 

If  Che  soil  is  homogeneous  and  can  be  claasified  as  clay  (cohesive) 
or  sand  (cohesionless ),  estimates  of  Eg  can  be  computed  using  values 
developed  by  Tersaghi  (11),  Brews  (9,  10),  and  others  based  on  experiments 
and  theoretical  relations.  The  range  of  values  provided  by  these  authors 
(see  next  sections)  aust  be  reduced  for  cyclic  and  group  effects  on  piles. 

It  should  be  noted  that  these  values  are  based  on  the  assuaption  that  the 
soil  has  linear  elastic  properties.  The  values  are  only  valid  for  siaple 
soil  conditions  and  aust  be  used  with  caution. 

Cohesive  Soils 

For  cohesive  soils  is  assumed  to  be  constant  with  depth, 

i.  -  v  a) 

“  Ej  (lft/»  (2) 

The  value  of  can  be  estimated  by  using  the  relation, 

Kj  -  aiSOq^  (3) 

where  a  is  a  paraaeter  on  a  1  foot  strip  that  varies  froa  0.32  to  0.52 
(Reference  10,  20),  generally  use  a  ■  0.4  (ft  *). 

2 

is  the  unconfined  compression  strength  of  clay  in  lbs/ft  . 

Thus, 

Kj  -  0.4(80^)  -  32q^  (4) 


Therefore,  ■  32 


(5) 


3  2 

Note  that  the  units  of  are  in  lbs/ft  with  ^  being  in  lbs/ft  and 


(6) 


B.  -  32qu  (6) 

Mote  that  the  units  of  Ea  in  the  above  equation  are  the  as 


as  those  of 


Reduction  Factors 

The  values  obtained  for  E#  wist  be  reduced  to  account  for  the  effects  of 
cyclic  loading  (Rj)  and  group  action  (Rj).  Thus, 

Use  Rj  *  1  for  initial  loading  and  ■  0.3  for  cyclic  loading  (References  9, 
10).  The  value  of  R2  can  be  obtained  froai  Table  D-l  (Reference  1)  given 
below. 


TABLE  0-1 


Values  of  Group  Effect  Factor  (R^) 


Pile  gpacing  in 
Direction  of  Loading 


Cohesionleass  8oils 

In  cohesionless  soils,  due  to  confineaent  effects,  Ba  inceases  with  depth. 
It  is  generally  assuaed  that  this  variation  is  linear. 

E#  -  n^Z  (8) 

when  n^  is  the  constant  of  horisontal  subgrade  reaction  and  Z  is  the  depth 
below  equivalent  ground  surface. 


Tersaghi  (Reference  11)  provided  the  following  empirical  relation  for 
obtaining  n^i 


1.3S 

where  A  i»  e  value  obtained  ea  a  function  of  the  relative  density  of  eoil  aa 
shown  in  Table  D-2  and  t  is  the  effective  unit  weight  of  sand.  The  units  of 
A  will  be  the  sane  as  those  of  f . 


TABLE  D-2 


Values  of  A  as  a  Function  of  the  Relative  Density  of  Send 
Relative  Density  H*  Value  of  A 

Loose  4-10  100-300 

Medium  10  -  30  300  *  1000 

Dense  30-50  1000  -  2000 


*H  is  the  number  of  blows  of  the  drop  weight  required  to  drive  the  sampling 
spoon  into  the  soil  for  a  distance  of  one  foot.  A  weight  of  140  lbs  end  a 
height  of  fall  of  30  inches  are  considered  standard. 


The  value  /  is  the  effective  unit  weight  of  the  sand.  If  the  piles  are  below 
the  water  table,  the  submerged  unit  weight  of  the  sand  should  be  used  in 
computing  the  value  of  u)x.  If  the  piles  extend  above  the  below  the  water 
table,  an  equivalent  height  of  submerged  sand  should  be  developed  above  the 
water  table.  The  depth  Z  then  must  be  measured  from  the  top  of  the 
equivalent  ground  surface.  Thus, 


watertable 


Seduction  factors  for  cycle  loading  (Rj )  and  group  loading  (S2)  are 
applied  in  the  ine  manner  for  cohesive  soils. 

Layered  or  Heterogeneous  Soil a 

In  general,  for  layered  or  heterogeneous  soil  deposits,  can  vary 
arbitrarily  with  depth  and  the  variation  is  difficult  to  be  represented  in  a 
mathematical  form.  For  these  cases,  p-y  curves  need  to  be  developed  along 
the  length  of  the  pile  and  computer  programs  as  explained  in  Appendix  C, 
Section  II,  C.2  can  be  used  to  compute  the  pile  stiffness  coefficients.  The 
p-y  curves  can  be  determined  based  on  laboratory  triaxial  testa  or  by 
measurements  of  the  behavior  of  instrumented  piles.  References  8,  20  and  23 
explain  the  developswnt  of  p-y  curves. 


For  homogeneous  layered  soils  that  are  cohesive  or  cohesionless,  Che 
appropriate  relations  given  in  the  preceding  sections  could  be  used.  For 
cohesive  soils,  the  Eg  values  are  only  dependent  on  the  unconfined 
compressive  strength  of  each  of  the  layers.  For  cohesionless  soils,  the  Eg 
values  are  dependent  on  n^  for  the  layer  and  the  effective  depth 
The  effective  depth  is  calculated  for  each  layer  as  shown: 


Z2 


X  2 

>7 


-  ( izi 
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(lz  1  ♦  <?2Z2 


Summary 

Typical  values  for  Eg  are  given  for  homogeneous  cohesive  and  cohesionless 
soils.  It  must  be  emphasized  that  the  pile  lateral  stiffness  coefficients 
can  be  directly  computed  only  when  Eg  is  a  constant  (such  as  for  a 
homogeneous  cohesive  soil)  or  linearly  varying  (as  for  homogeneous 
cohesive less  soil).  In  all  other  cases  including  layered  soils  and  soils 
that  cannot  be  classified  as  cohesive  or  cohesionless,  the  pile  stiffness 
coefficients  can  be  calculated  by  using  computer  programs  such  as  those 
elaborated  in  Appendix  C,  Section  II,  C.2  of  this  report. 
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APPENDIX  E 
EXAMPLE  PROBLEMS 


EXAMPLE  PROBLEMS  1  THROUGH  5 

Two-dimensional  Problems 

Hrennikoff’s  Example 

Comparison  of  Computer  Output  with 
the  Elastic  Center  Method 

Hrennikoff's  calculations  and  a 
Common  Analytical  Method 

Example  problems  1  through  5  illustrate  how  the  computer  can  be  used  to  analyze 
two-dimensional  pile  foundation  problems.  The  following  examples  were  taken 
from  Hrennikoff's  paper  entitled  "Analysis  of  Pile  Foundations  with  Batter  Piles", 
published  in  the  Transactions  of  the  American  Society  of  Civil  Engineers,  Vol. 

76,  No.  1,  Paper  No.  2401,  Jan.  1950,  pp.  123-126. 

The  computer  results  are  compared  with  Hrennikoff's  results  as  well  as  with 
other  hand  computation  methods  commonly  used  by  civil  engineers. 

The  physical  pile  layout  for  example  problems  1  through  5  is  shown  in  Figure  El. 


Example  Problem  No.  1  compares  the  results  obtained  by  the  Computer  Method  with 
those  obtained  by  the  Elastic  Center  Method  assuming  the  soil  offers  no  lateral 
support;  in  other  words,  the  subgrade  modulus  is  zero.  The  Computer  results 
agree  closely  with  the  Elastic  Center  Method  results.  A  description  of  the 
Elastic  Center  Method  can  be  found  in  "Substructure  Analysis  and  Design"  by 
Dr.  Paul  Andersen.  This  method  is  limited,  however,  to  pile  groups  consisting 
of  hinged  piles  arranged  in  two  subgroups  whose  centerlines  intersect. 


El 


Example  Problem  No.  2  compares  the  results  obtained  by  the  Computer  Method 
with  the  results  obtained  by  Manual  Calculations  as  presented  in  Hrennikoff's 
paper  for  case  2a  (very  weak  soil  with  a  subgrade  modulus  of  3.123  psi) .  The 
Computer  results  agree  closely  with  Hrennikoff's  calculations. 

Example  Problem  No.  3  compares  the  results  obtained  by  the  Computer  Method  with 
the  results  obtained  by  Manual  Calculations  as  presented  in  Hrennikoff's  paper 
for  Case  4a  (weak  soil  with  a  subgrade  modulus  of  31.23  psi).  The  Computer 
results  agree  closely  with  Hrennikoff's  calculations. 

Example  Problem  No.  4  compares  the  results  obtained  by  the  Computer  Method  with 
the  results  obtained  by  Manual  Calculations  as  presented  in  Hrennikoff's  paper 
for  Case  6a  (medium  soil  with  a  subgrade  modulus  of  312.3  psi).  Again,  the 
Computer  results  agree  closely  with  Hrennikoff's  calculations. 

Example  Problem  No.  5  compares  the  results  obtained  by  the  Computer  Method  with 
the  results  obtained  by  a  common  Analytical  Method  for  two  different  load 
conditions.  A  description  of  the  common  Analytical  Method  can  be  found  in 
"Foundation  Engineering"  by  Ralph  Peck,  Walter  E.  Hanson,  and  Thomas  H. 

Thornbum,  and  in  "Foundation  Design"  by  Wayne  Teng.  Example  Problem  No.  5 
demonstrates  that  the  individual  pile  forces  obtained  by  the  common  Analytical 
Method  are  approximate  and  may  or  may  not  agree  closely  with  the  results  obtained 
by  the  Computer  Method.  A  subgrade  modulus  of  312  psi  was  used  for  this  example. 
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Example  Problem  1 


1 


Two-dimensional  problem 
Hrennikoff's  example 
no  lateral  soil  support 


Properties 

E  =  0.15  x  107  psi 

Degree  of  fixity  =  0.0 

Mod  .  =  0.1  »0 
sub 

Pile  resistance  =  0.5 

I  =  322.06  in.4 

Participation  factor 

for  torsion  =  0.0 

I  =  322.06  in.4 
y 

Area  =  63.5  in.2 

Length  =  30  ft 

Torsion  modulus  =  0.0 

Loading 

Case 

Q3 

(kips) 

Q5 

(kip-ft) 

1 

-39.375 

113.1 

173.4 

Properties  and  loading  conditions 
for  example  problem  1 


1  example  problem  NO  1 

2  HKENNIKOFF  EXAMPLE 

3  2 

4  5  11 

5  1  .1 

4  1  5  30.000  3 

7  9.000 

8  4 

9  1500000.000 

10  2 

11  0.  1.0  0.  0. 

12  1 

13  82.000  40.000 

14  2 

15  1  3  -3.000 
14  4  5  0. 

17  -5.000  -2.5  0.  3.000  7.000 

18  -39.375  113.1  173.4 
E>FILE 

STR0M5  EDITED. 

OREPLACE.STROH5 
OOLDr  CORPS  A1N-CECELB 
DCALL.C0RPS.X0034 

JOB  WAITING. 


INPUT  DATA  FILE  NAME  IN  7  CHARACTERS  OR  LESS.  HIT  A 
CARRIAGE  RETURN  IF  INPUT  DATA  MILL  COME  FROH  TERMINAL. 


DSTR0H5 


THIS  PROGRAM  GENERATES  THE  FOLLOWING  TABLES! 


TABLE  NO. 
1 
2 

3 

4 

5 
4 

7 

8 


CONTENTS 

PILE  AND  SOIL  DATA 

PILE  COORDINATES  AND  BATTER 

STIFFNESS  AND  FLEXIBILITY  MATRICES  FOR  THE 

STRUCTURE  AND  COORDINATES  OF  ELASTIC  CENTER 

APPLIED  LOADS 

STRUCTURE  DEFLECTIONS 

PILE  DEFLECTIONS  ALONG  PILE  AXIS 

PILE  FORCES  ALONG  PILE  AXIS 

PILE  FORCES  ALONG  STRUCTURE  AXIS 


E5 


INPUT  THE  NUMBERS  OF  THE  TABLES  FOR  WHICH  YOU  WANT  THE  OUTPUT, 
SEPARATE  THE  NUMBERS  WITH  COMMAS. 

I>li2i3i4iStir7i8 


INPUT  A  FILENAME  FOR  TABLE  8  IN  7  CHARACTERS  OR  LESS 
IF  TOU  WANT  TO  USE  THIS  INFORMATION  FOR  A  NEW  RUN 
HIT  A  CARRIAGE  RETURN  IF  YOU  DO  NOT  WANT  THIS  FILE. 


I) 


INPUT  A  FILE  NAME  FOR  OUTPUT  IN  7  CHARACTERS  OR  LESS. 

HIT  A  CARRIAGE  RETURN  IF  OUTPUT  IS  TO  BE  PRINTED  ON  TERMINAL. 


I> 


EXAMPLE  PROBLEM  NO  1 
HRENNIKOFF  EXAMPLE 


NO.  OF  PILE  ROWS  -  5  B  MATRIX  IS  CALCULATED  FOR  EACH  ROW 


tmmmmmmmuumtmtmmmmummmmmmmm 


1.  TABLE  OF  PILE  AND  SOIL  DATA 


PILE  NUMBERS 


E  =  . 

I5EF07  PSI  IX 

=  322.04  IN«4  IT  = 

322.04  IN«4 

AREA  = 

43.4  INM2 

X  *  9.00  IN  Y  = 

9.00  IN 

LENGTH 

=  30.0  FEET 

ES  =  .100 

K1  = 

.4107  K2  * 

1.0000  K3  =  0.0000 

K4  = 

0.0000  KS  = 

0.0000  K4  =  0.0000 

LENGTH  OF  PILES  (  30.00  FEET)  IS  INSUFFICIENT 
FOR  PILE  GROUP  -  1  MINIMUM  ACCEPTABLE  LENGTH  IS  B7.B8  FEET 

FOR  SEMI-INFINITE  BEAM  ON  ELASTIC  FOUNDATION 


ALLOWABLES:  COMPRESSIVE  LOAD  :  82.000  KIPS 

TENSILE  LOAD  =  40.000  KIPS 


THE  B  MATRIX  FOR  PILES  1  THROUGH  5  IS 


•184E402  0.  0, 
0.  .245E+04  0. 
0.  0.  0. 


ttmmummmmmmmtmtmmuummttmtmmtmmt 


2.  inilt  OF  PILE  COORDINATES  AND  BATIER 


PILE  ROW 

BATTER 

U1  <FT) 

1 

-3.00 

-5,000 

2 

-3.00 

-2.500 

3 

-3.00 

0.000 

4 

VERTICAL 

3,000 

5 

VERTICAL 

7.000 

mummttmummtmmummummtmummtmmmm 
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3.  STIFFNESS  MATRIX  S  FOR  THE  STRUCTURE 


.794EF05  — .239E+06  -.714E+07 
-.239EM  .125EW  -,103£m 
-.71AET07  -.103EM8  .329E+10 


3A  FLEXIBILITY  MATRIX  F  FOR  THE  STRUCTURE 


.127E-03  .274E-04  .343E-04 
.274E-04  .471E-05  .B07E-07 
.343E-04  .807E-07  .135E-0B 


COORDINATES  OF  ELASTIC  CENTER 
EC1  =  .003  EC2  =  -.002 


tmtm  LOADING  CONDITION  1  Mttttt 


4.  MATRIX  OF  APPLIED  LOADS  0  <KIPS  I  FEET) 
01  03  05 

-39.375  113,106  173.400 


mmmuumnmmtttmmmmttmutmummmummm 


E8 


it,  ■- 


-.116E+01  150E+00  -.234E-02 


mmnuutmmmmmmmtmmmmtmmummtmmm 


6.  PILE  DEFLECTIONS  ALONG  PILE  AXIS  (INCHES) 


PILE 

XI 

X  3 

XS 

1 

-.119E+01 

.890E-01 

-.236E-02 

2 

-.117E+01 

.156EF00 

-.236E-02 

3 

-.115E+01 

.223E+00 

-.236E-02 

4 

-.116E+01 

-.455E-01 

-.236E-02 

5 

-.116E+01 

.478E-01 

-.236E-02 

ttuumumtmmmmmmmmtmttmmtmmmmmttm 


7.  PILE  FORCES  ALONG 

PILE  AXIS  (KIPS  1  FT) 

PILE 

FI 

F3 

F5 

FAILURE 

BU  CO  TE 

1 

-.022 

23.589 

0.000 

F 

2 

-.022 

41.394 

0.000 

F 

3 

-.021 

59.199 

0.000 

F 

4 

-.022 

-17.359 

0.000 

S 

-.022 

12.670 

0.000 

TOTAL 

NO.  FAILURES  = 

3 

LOAD  CASE  1 

mmmttmmmtmummtmmmtmmmmummmmm 


8.  PILE  FORCES  ALONG  STRUCTURE  AXIS  (KIPS  i  FEET) 


PILE 

FI 

F3 

F5 

1 

-7.481 

22.371 

0.000 

2 

-13.111 

39.243 

0.000 

3 

-18.741 

56.155 

0.000 

4 

-.022 

-17.359 

0.000 

5 

-.022 

12.470 

0.000 

Force  Diagram 
Forces  resolved m 
d I  rechon  of  pile  groups 


Consider  Piles  /,  Zd  o  are  pile  group  A 
and  Pdes  4-$ 5  are  pile  group  S 

Thm  by  the  e io s he  center  'method, 

'°A  -  4-  r  M 

nA  8rz 

and  Pe  =.  ft 6  +J2 _  Pi 

na  8rz 

cohere  '■  Ra?Rb  are  force  components  m  directions  of 
A  ?  S  pile  groups 

Da  i'Hd  ore  number  of  piles  j/)  A?  &  pile  QroupS 

r  is  the  distance  Prom  the  dnz  he  center  to the pile 

M  15  the  moment  of  the  applied  loud  about 

the  ehshc  center. 


£rz  =  2  (zdsD1  *  z(z.o)2'  -  13.05 

Pile  i  =  JZ ±_  _  /.  E  6/9-76.)  2. 35  _  23.<*'K 

3  13.05 

Pile  Z  =  /?-.!_  _  4/.  33  < 

3 


Pile  3  =  iZ 4  t  /.  2.  ( H3.IC0)  2-35 

3  /  9. 05 

P< /S  4  -  -  -3  _  _/. 2 / //9-7g)E 

a  ""‘/9.0  s* 

Pile  5  =  -  -f  /.2  (//9.7<o)Z 
Z  ’3.0S 


^  53.0<ok 
=  -fV-59  < 

=  /2-59< 


Comparison  of  results  obtained  by  the  Compuhe 
r> nerhoa'  w/ih  the  results  obtomed 6y  the 
elastic  center  method. 


P/Je  Forces  atony  Pile  Axis 


Pile 

do. 


<  ofiipnter 
Output 

Ps 

Oops) 

Eloshc 

Center 

Fs 

C&p  s ) 

£3.S&9 

23.0! 

41.  394- 

41.  33 

5  9-/99 

59.06 

-17.353 

-  n.  59 

IZ.07O 

/  2.59 

e: 


Example  Problem  2 

Two-dimensional  problem, 
Hrennikoff 's  example 
case  2a  (very  weak  soil) 


Properties 

E  =  0.15  *  10^  psi 

Degree  of  fixity  = 

0.0 

Modsub  =  3-123  Psi 

Pile  resistance  = 

0.5 

Ix  =  322.06  in.U 

Participation  factor 

for  torsion  = 

0.0 

I  =  322.06  in.14 

Torsion  modulus  = 

0.0 

Area  =  63.5  in.2 

Length  =  30  ft 

Loading 

Case 

Q1 

(kips) 

Q5 

(kip-ft) 

1 

-39.375 

113.1 

173. 

Properties  and  loading  conditions  for 
example  problem  2 
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1  EXAMPLE  PROBLEM  NO  2 

2  HRENNIKOFF  EXAMPLE 

3  2 

AS11 

5  1  3.123 

6  1  S  30.000  3 

7  9.000 

8  4 

9  15 00000.000 

10  2 

11  0.  .5  0.  0. 

12  1 

13  82.000  40.000 

14  2 

15  1  3  -3.000 

16  4  5  0. 

17  -5.000  -2.5  0.  3.000  7.000 

18  -39.375  113.1  173.4 
E>FILE 

STR0M5  EDITED. 

DREPLACE.STR0M5 
DOLD.  CORPS. 

ERROR  IN  ARGUMENT. 

DOLD.  C0RPS/UN-CECEL8 
OCALL  .CORPS.  X0034 


INPUT  DATA  FILE  NAME  IN  7  CHARACTERS  OR  LESS,  HIT  A 
CARRIAGE  RETURN  IF  INPUT  DATA  HILL  CONE  FROM  TERMINAL. 


I>STR0N5 


THIS  PROGRAM  GENERATES  THE  FOLLOHING  TABLES! 


TABLE  NO.  CONTENTS 

1  PILE  AND  SOIL  DATA 

2  PILE  COORDINATES  AND  BATTER 

3  STIFFNESS  AND  FLEXIBILITY  MATRICES  FOR  THE 
STRUCTURE  AND  COORDINATES  OF  ELASTIC  CENTER 

4  APPLIED  LOADS 

5  STRUCTURE  DEFLECTIONS 

6  PILE  DEFLECTIONS  ALONG  PILE  AXIS 

7  PILE  FORCES  ALONG  PILE  AXIS 

8  PILE  FORCES  ALONG  STRUCTURE  AXIS 
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INPUT  THE  NUMBERS  OF  THE  TABLES  FOR  WHICH  TOU  WANT  THE  OUTPUT. 
SEPARATE  THE  NUMBERS  WITH  COMMAS. 

I>l»2i3.-4i5.Bf7»8 


INPUT  A  FILENAME  FOR  TABLE  8  IN  7  CHARACTERS  OR  LESS 
IF  YOU  WANT  TO  USE  THIS  INFORMATION  FOR  A  NEW  RUN 
HIT  A  CARRIAGE  RETURN  IF  YOU  DO  NOT  WANT  THIS  FILE. 


tAAflrlE  PROBLEM  NO  2 
HRENNIKOFF  EXAMPLE 

NO.  OF  PILE  ROWS  =  5  B  MATRIX  IS  CALCULATED  FOR  EACH  ROW 

tmmmMmmmmmmmttmmmmmwmmmmmttt 


( 
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1.  TABLE  OF  PILE  AND  SOIL  DATA 


PILE  NUMBERS 


1  5  £  1  .15EW7  PSI  IX  =  322.04  INttA  IT  =  322.04  INM4 

AREA  =  43.4  IN»2  X  =  9.00  IN  T  =  9.00  IN 

LENGTH  =  30.0  FEET  ES  =  3.123 

K1  =  .4107  K2  =  .5000  K3  =  0.0000 

K4  =  0.0000  K5  =  0.0000  K4  =  0.0000 


LENGTH  OF  PILES  (  30.00  FEET)  IS  INSUFFICIENT 
FOR  PILE  GROUP  -  1  HINIHUH  ACCEPTABLE  LENGTH  IS  37.17  FEET 

FOR  SEMI-INFINITE  BEAM  ON  ELASTIC  FOUNDATION 


ALLOWABLES!  COMPRESSIVE  LOAD  =  82.000  KIPS 

TENSILE  LOAD  =  40.000  KIPS 


THE  B  MATRIX  FOR  PILES  1  THROUGH  5  IS 


.244E+03  0.  0. 
0.  .133EF04  0. 
0.  0.  0. 


tmmmutmmmmmmmmmmtmmtmmmmmmm 


2.  TABLE  OF  PILE  COORDINATES  AND  BATTER 


PILE  ROW 

BATTER 

U1  (FT) 

1 

-3,00 

-5.000 

2 

-3.00 

-2.500 

3 

-3.00 

0.000 

4 

VERTICAL 

J.000 

5 

VERTICAL 

7.000 

E17 


3.  STIFFNESS  MATRIX  S  FOR  THE  STRUCTURE 


•409E+05  -.119E406  -.357E407 
-.119E+06  .623E+04  -.517E+07 
-.357E407  -.517E407  .1HE410 


3A  FLEXIBILITY  MATRIX  F  FOR  THE  STRUCTURE 


. 193E-03  .414E-04  .550E-04 
.414E-04  .105E-04  .123E-0A 
.550E-06  .123E-04  .219E-08 


COORDINATES  OF  ELASTIC  CENTER 
EC1  =  .003  EC2  *  -.002 


mum  L0ADIN6  CONDITION  1  mttttt 


tmmmtmmmmmmmmnmmtmmumummmmm 


A.  MATRIX  OF  APPLIED  LOADS  Q  (KIPS  t  FEET) 
01  03  05 

-39.375  113.100  173.400 


mtmtkmmutMmmtmummmummnmmmmmttm 


S.  STRUCTURE  DEFLECTIONS  (INCHES) 


D1  D3  D5 

-.177E+OI  - >183E400  -.315E-02 


tmtmnmtmtmmumtmtutmttummtmttntmmttmu 


t.  PILE  DEFLECTIONS  ALONG  PILE  AXIS  (INCHES) 
PILE  XI  X3  X5 

1  -.180EA01  .207E+00  -.315E-02 

2  -.177EA01  .290E+00  -.315E-02 

3  -.17AE401  .384E400  -.315E-02 

A  -.177E+01  -.492E-01  -.315E-02 

5  -.177EA01  .821E-01  -.315E-02 


mwmmmtmmtmmmmmtutmtmmmmmmmtmt 


7,  PILE  FORCES  ALONG  PILE  AXIS  (KIPS  1  FT) 


PILE 

FI 

F3 

F  5 

FAILURE 
BU  CO  TE 

1 

->AA2 

27.395 

0.000 

F 

2 

-.A35 

39.282 

0.000 

F 

3 

->A27 

51.170 

0.000 

F 

A 

-.A3 6 

-P.107 

0.000 

5 

->A3  6 

10.881 

0.000 

TOTAL  NO,  FAILURES  =  3  LOAD  CASE  1 


umutmmmmtmmmmtmmmtmnmmtmmmttunt 
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B.  PILE  FORCES  ALONG  STRUCTURE  AXIS  (KIPS  I  FEET) 


PILE 

FI 

F3 

F5 

1 

-9.082 

25.849 

0.000 

2 

-12.835 

37.129 

0.000 

3 

-16.587 

43.403 

0.000 

4 

-.436 

-9.147 

0.000 

5 

-.434 

10.381 

0.000 

SUH 

-39.375 

113.1*0 

173.400 

ummmmmmtmmtmmmmtmttmmmttnmmmtm 
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Commit  er 

Out out 

Hrennikoff '  s 

Examole 

Pile 

F1 

F3 

F1 

F3 

No. 

(kios) 

(kios ) 

(kios) 

(kips) 

1 

0 .  UU2 

27.395 

0.1*1* 

27.5 

2 

0.1*35 

39.282 

0.1*3 

39-3 

3 

0.1*27 

51.170 

0.1*3 

51.0 

1* 

0.1*36 

-9-167 

0.1*3 

-9.0 

5 

0.1*36 

10.881 

0.1*3 

10.9 

E21 


Example  Problem  3 


Two-dimensional  problem, 
Hrennikoff's  example  • 
case  ka  (weak  soil) 


Properties 


E  =  0.15  x  107  psi 

Modsub  =  31-230  psi 
I  =  322.06  in. 

X 

I  =  322.06  in.u 

y  o 

Area  =  63-5  in. 

Length  =  30  ft 


Degree  of  fixity  =0.0 

Pile  resistance  =1.0 

Participation  factor 

for  torsion  =0.0 

Torsion  modulus  =0.0 


Loading 

Case 


"‘I 

(kips) 


-39-375 


s 

(kips) 


(kij>-ft) 


113.1 


173.  U 


Properties  and  loading  conditions  for 
example  proolem  3 


1  EXrtfr'Lt  PKUBltn  ftU  3 

2  HRENNIKOFF  EXAMPLE 

3  2 

4  5  11 

5  1  31.23 

4  1  5  30.000  3 

7  9.000 

8  4 

9  1500000.000 

10  2 

11  0.  1.0  0.  0. 

12  1 

13  82.000  40.000 

14  2 

15  1  3  -3. COO 
14  4  5  0. 

17  -5.000  -2.5  0.  3.000  7.000 

18  -39.375  113.1  173.4 
E>FILE 

STROMS  EDITED. 

OREPLACE.STRQM5 
DOLDiCORPCS 
ERROR  IN  ARGUMENT. 
C>OLD»CORPS/UN=C£CELB 
DCALLrC0RPSrX0034 


INPUT  DATA  FILE  NAME  IN  7  CHARACTERS  OR  LESS.  HIT  A 
CARRIAGE  RETURN  IF  INPUT  DATA  HILL  COKE  FROM  TERMINAL. 


DSTR0M5 


THIS  PROGRAM  GENERATES  THE  FOLLOWING  TABLES! 

TABLE  NO.  CONTENTS 

1  PILE  AND  SOIL  DATA 

2  PILE  COORDINATES  AND  BATTER 

3  STIFFNESS  AND  FLEXIBILITY  MATRICES  FOR  THE 
STRUCTURE  AND  COORDINATES  OF  ELASTIC  CENTER 

4  APPLIED  LOADS 

5  STRUCTURE  DEFLECTIONS 

4  PILE  DEFLECTIONS  ALONG  PILE  AXIS 

7  PILE  FORCES  ALONG  PILE  AXIS 

6  PILE  FORCES  ALONG  STRUCTURE  AXIS 


INPUT  THE  NUMBERS  OF  THE  TABLES  FOR  WHICH  YOU  WANT  THE  OUTPUT. 
SEPARATE  THE  NUMBERS  WITH  COMMAS. 

I>1»2»3M.5.4»7»8 


INPUT  A  FILENAME  FOR  TABLE  8  IN  7  CHARACTERS  OR  LESS 
IF  YOU  WANT  TO  USE  THIS  INFORMATION  FOR  A  NEW  RUN 
HIT  A  CARRIAGE  RETURN  IF  YOU  DO  NOT  WANT  THIS  FILE. 


I> 


EXAMPLE  PROBLEM  NO  3 
HRENNIKOFF  EXAMFLE 


NO.  OF  PILE  ROWS  *  5  B  MATRIX  IS  CALCULATED  FOR  EACH  ROW 


mmmmmmtmtmmtuutmmmnmtmmmtummmt 


1.  TABLE  Of  PILE  AND  SOIL  BATA 


p;le  numbers 


i 


5  E  =  .15E407  PSI  IX 
AREA  =  A3, 6  IN**: 

LENGTH  =  30.0  FEET 

K1  =  .4107  K2  = 

K4  =  0.0000  K5  = 


=  322.06  IN**4  IY  = 

X  =  9.00  IN  Y  = 

ES  =  31.230 
1.0000  K3  =  0.0000 

0,0000  KA  =  0.0000 


322.06  IN**4 
9.00  IN 


allowables:  COHPRESSM  LOAD  =  S2.000  KIPS 

TENSILE  LOAD  =  40.000  KIPS 


THE  B  MATRIX  FOR  PILES  1  THROUGH  5  IS 


.138E404  0.  0. 
0.  .2A5E406  0. 
0.  0.  0. 


ttmmnmtmnmmuttmmmmmmmtmmmmmmm 


2.  TABLE  OF  PILE  COORDINATES  AND  BATTER 


PILE  RON 

BATTER 

U1  <FT) 

1 

-3.00 

-5,000 

2 

-3.00 

-2.500 

3 

-3.00 

0.000 

4 

VERTICAL 

3.000 

S 

VERTICAL 

7. MO 
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3.  STIFFNESS  MATRIX  S  FOR  THE  STRUCTURE 


.84QEF05  -.237EM  -.712E+07 
-.237EE06  . 125E+07  -.1C3E108 
-.712EF07  I03E+08  .329ET10 


3A  FLEXIBILITY  MATRIX  F  FOR  THE  STRUCTURE 


.664E-04  . 142E-04  .183E-0* 
• 142E-04  .3S3E-05  .429E-07 
.188E-0S  .429E-07  .847E-09 


COORDINATES  OF  ELASTIC  CENTER 
EC1  --  .003  EC2  -  -.002 


mum  LOADING  CONDITION  i  ttmtu 


mmnmmmmmmmmmmmmmmtmmmmmmm 

4.  MATRIX  OF  APPLIED  LOADS  0  (KIPS  I  FEET) 

Q1  03  05 

-39.375  113.100  173.400 

mmumtmmmtmmttmmmtmmmmmummmmm 
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5.  STRUCTURE  DEFLECTIONS  (INCHES) 


D1  03  D5 

-♦616E+00  -.332E-01  -.805E-03 


4.  PILE  DEFLECTIONS  ALONG  PILE  AXIS  (INCHES) 
FILE  XI  X3  XS 


I  -.410E+00 

.117EFOO 

-.805E-03 

2  -.402E400 

.140E400 

-.805E-03 

3  -.595E+00 

.163E400 

-.805E-03 

4  -.414E400 

-.421E-02 

-.805E-03 

5  -.414E+00 

•344E-01 

-.805E-03 

mmmmmummmtmtmummmmtmttmmmummt 


7.  PILE  FORCES  ALONG  PILE  AXIS  (KIPS  t  FT) 


PILE 

FI 

F3 

F5 

FAILURE 
BU  CO  TE 

1 

-.845 

31.132 

0.000 

F 

2 

-.834 

37.2C4 

0.000 

F 

3 

-.824 

43.273 

0.000 

F 

4 

-.853 

-1.117 

0.000 

S 

-.853 

9.124 

0.000 

TOTAL  NO.  FAILURES  =  3  LOAD  CASE  1 


mmmtmmmtmmmmmmmtmumtmttmtmmmm 
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8,  PILE  FORCES  ALONG  STRUCTURE  AXIS  (KIPS  t  FEET) 


PILE 

FI 

F3 

F5 

1 

-10.646 

29.267 

0.000 

2 

-12.553 

35.031 

0.000 

3 

-14.467 

40.795 

0.000 

•4 

-.853 

-1.117 

0.000 

5 

-.853 

9.124 

0.000 

SUN 

-39.375 

113.100 

173.400 
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Hrennikoff *s 


Computer 

Outuut 

Examule 

Pile 

No . 

F1 

(kins ) 

F3 

(kins) 

F1 

(kios) 

F3 

(kios ) 

1 

O.8U5 

31-132 

0.8U 

31.2 

2 

0.83U 

37.201+ 

0.83 

37.2 

3 

0.82U 

2+3.276 

0.82 

1+3.2 

k 

0.853 

-1.117 

0.85 

-1.0 

5 

0.853 

9.12U 

0.85 

9.1 
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Example  Problem  4 


Two-dimensional  problem, 
Hrennikoff's  example 
case  6a  (medium  soil) 


Properties 

E  =  0.15  *  10^  psi 

Degree  of  fixity  = 

0.0 

Mod  .  =  312.30  psi 
sub 

Pile  resistance  = 

1.00 

I  =  322.06  in. 

Participation  factor 

J. 

for  torsion  = 

0.0 

I  =  322.06  in. 

y  2 

Area  =  63.5  in. 

Torsion  modulus  = 

0.0 

Length  *  30  ft 

Loading 

Case 

Q1 

(kips) 

Q5 

(kip-ft) 

1 

-39-375 

113.1 

173.  »* 

Properties  and  loading  condition  for 
example  problem  4 


i 

s 


2  HRENN1MFF  EXAMPLE 

3  2 

4  5  11 

5  1  312.3 

6  1  5  30.000  3 

7  9.000 
6  4 

9  1500000.000 

10  2 

11  0.  1.0  0.  0. 

12  1 

13  82.000  40.000 

14  2 

15  1  J  -3.000 
14  4  5  0. 

17  -5.000  -2.5  0.  3.000  7.000 

18  -39.375  113.1  173.4 
f>FIL£ 

STR0M5  EDITEO. 

DKEPL  ACE.  STROMS 
C>OL  D . CORPS/UN=CEC£L  6 
DCALL.  CORPS  .X0034 


INPUT  DATA  FILE  NAME  IN  7  CHARACTERS  OR  LESS.  HIT  A 
CARRIAGE  RETURN  IF  INPUT  DATA  UILL  COME  FROM  TERMINAL. 


DSTR0M5 


THIS  PROGRAM  GENERATES  THE  FOLLOHING  TABLES.' 

TABLE  NO.  CONTENTS 

1  PILE  AND  SOIL  DATA 

2  PILE  COORDINATES  AND  BATTER 

3  STIFFNESS  AND  FLEXIBILITY  MATRICES  FOR  THE 
STRUCTURE  AND  COORDINATES  OF  ELASTIC  CENTER 

4  APPLIED  LOADS 

5  STRUCTURE  DEFLECTIONS 

6  PILE  DEFLECTIONS  ALONG  PILE  AXIS 

7  PILE  FORCES  ALONG  PILE  AXIS 

8  PILE  FORCES  ALONG  STRUCTURE  AXIS 
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INPUT  THE  NUMBERS  OF  THE  TABLES  FOR  WHICH  YOU  KANT  THE  OUTPUT. 
SEPARATE  THE  NUMBERS  KITH  COMMAS. 

I>l»2f3fA»5»6.7»B 


INPUT  A  FILENAME  FOR  TABLE  8  IN  7  CHARACTERS  OR  LESS 
IF  YOU  WANT  TO  USE  THIS  INFORMATION  FOR  A  NEW  RUN 
HIT  A  CARRIAGE  RETURN  IF  YOU  DO  NOT  WANT  THIS  FILE. 


I> 


EXAMPLE  PROBLEM  NO  4 
HRENNIKOFF  EXAMPLE 

NO.  OF  PILE  ROUS  -  5  B  MATRIX  IS  CALCULATED  FOR  EACH  RDU 


ttmmmmmmmumummmmmmmmmtmummm 
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1.  TABLE  OF  PILE  AND  SOIL  DATA 


PILE  NUMBERS 


1  5  £:  *  15E+07  PSI  IX  =  322,06  INW4  IY  =  322.06  IN«4 

AREA  =  63.6  IN«2  X  =  9.00  IN  Y  =  9.00  IN 

LENGTH  =  30.0  FEET  ES  =  312.300 

K1  =  .4107  K2  -  1.0000  K3  =  0.0000 

K4  =  0.0000  KS  =  0.0000  K6  =  0.0000 


ALLOWABLES:  COMPRESSIVE  LOAD  =  82.000  KIPS 

TENSILE  LOAD  =  40.000  KIPS 


THE  B  MATRIX  FOR  PILES  1  THROUGH  5  IS 


. 779E+04  0.  0. 
0.  .265E+06  0. 
0.  0.  0. 


tmmummmmtmmmmmmmmmtmtmmmmmm 


2.  TABLE  OF  PILE  COORDINATES  AND  BATTER 


PILE  ROU 

BATTER 

U1  (FT) 

1 

-3.00 

-5.000 

2 

-3.00 

-2.500 

3 

-3.00 

0.000 

4 

VERTICAL 

3,00" 

X 

IMKTlCkL 

7.  AM 
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3,  STIFFNESS  MATRIX  S  FOR  THE  STRUCTURE 


.116EF04  -.232E404  -.495E407 
-.232E40A  .125E+07  -.103EI0G 

&95E+07  - . 103E+08  .329E410 


3A  FLEXIBILITY  MATRIX  F  FOR  THE  STRUCTURE 


.205E-04  .427E-05  .S44E-07 
•427E-05  .171E-05  .144E-07 
.56AE-07  .144E-07  .448E-09 


COORDINATES  OF  ELASTIC  CENTER 
EC1  =  .003  EC2  =  -.002 


tUtUtt  LOADING  CONDITION  1  I* MM 


tmtmmtmmmuumtmMtuMmmttmtmmmmmtm 


4.  MATRIX  OF  APPLIED  LOADS  Q  (KIPS  1  FEET) 
01  Q3  OS 

-39.373  113.100  173.400 


mmmumutmmmnmtmtmMttmmmmmttmutmtt 


E34 


5,  STRUCTURE  DEFLECTIONS  (INCHES) 
D1  93  DS 

-.207E+00  .553E-01  .348E-03 


iiminittmuttiimmmumutiittiiiittmtmmiuutmmiti 


6.  FILE  DEFLECTIONS  ALONG  PILE  AXIS  (INCHES) 


XI 

X3 

X5 

-.172ET00 

.139EW0 

.348E-03 

-.175E400 

. 128E+00 

.348E-03 

-.179E+00 

•  U8EF00 

.348E-03 

-.207E+00 

.420E-Q1 

•368E-03 

-.207E+00 

.243E-01 

.368E-03 

mMnmntumunmmtmtmmmnmtmmmtmmmmt 


7.  FILE  FORCES  ALONG  PILE  AXIS  (KIPS  l  FT) 


PILE 

FI 

F3 

F5 

I 

-1.338 

36.790 

0.000 

2 

-1.345 

34.014 

0.000 

3 

-1.392 

31.237 

0.000 

A 

-1.411 

11.137 

0.000 

5 

-1.411 

4.454 

0.000 

FAILURE 
BU  CO  TE 


TOTAL  NO.  FAILURES  = 


LOAD  CASE  1 


mmutumtttmmtmmtmtttmtmmmnmtmmtmmm 
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8.  PILE  FORCES  ALONG  STRUCTURE  AXIS  (KIPS  1  FEET) 


PILE 

FI 

F3 

F5 

1 

-12.903 

34.479 

0.000 

2 

-12.051 

31.834 

0.000 

3 

-11.199 

29.194 

0.000 

A 

-1.411 

11.137 

0.000 

5 

-1.411 

4.454 

0.000 

SUN 

-39.375 

113.100 

173.400 

Results  and  calculations 

Manual  calculations  for  this  example  are  presented  in 
Hrennikoff's  paper,  case  6a.  The  computer  results  shown 
agree  closely  with  the  classical  method  results.  For  example,  a  com¬ 
parison  of  the  horizontal  forces  in  each  pile  is  shown  "below: 

F^  (kips)  from 


Pile 

Computer 

Hrennikoff's 

No. 

Out  put 

Examnle 

1 

1.338 

1.3** 

2 

1.365 

1.37 

3 

1.392 

1.39 

1* 

1.611 

1. 6l 

5 

1.611 

l.6l 

The  vertical  pile  forces  also  agree  closely  and  are  shown 
"below: 


F^  (kips)  from 


Pile 

Computer 

Hrennikoff’s 

No. 

Outout 

Examole 

1 

36.790 

36.8 

2 

3U.01U 

3U.0 

3 

31.237 

31.2 

k 

11.137 

11.1 

5 

6.45U 

6.5 
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IzKOmp/e  Problems  /-4 

P/(e  Force  ~  F5  (tips)  A  loro  R/eAxis 
VS  Sob  groc/e  Mo c/o /os  Cpst) 


Pi  le 
hjo 


I  Subgrode  Modulus  Voloe  ( 

psi) 

0 

3.123 

3/. 23 

3/2.3 

^3.sa9K 

27.395* 

3/./32e 

36.790* 

4/.  394* 

39.262* 

37-204* 

34.0/4* 

5  9./ 99" 

5i.no * 

4-3.27  cP 

3/.237* 

-/7.3S9* 

-9. 16  7* 

-t.m*- 

//.  /37* 

/2.670* 

/0. 68/* 

9./24* 

6.454* 

Pile  Force -F3  (tips') 


Example  Problem  5 
Two-dimensional  Problem 

Batter  and  vertical  piles  supporting  a  wall  foundation  with  constant 
soil  moduli 


Properties: 

E  *  .15  x  10^  psi  (wood) 

Mod  .  -  312 
sub  , 

I  *  322  in 

x  4 

I  -  322  in 

y  2 
Area  “  63.5  in 

Length  “  30  ft 

Degree  of  fixity  “0.0 

Pile  resistance  “1.0 

Participation  factor  for  torsion  “0.0 

Torsion  modulus  »  0.0 

Loading: 


Loading 

Case 


THIS  PROGRAM  GENERATES  THE  FOLLOWING  TABLES 


TABTE  NO.  CONTENTS 

1  PILE  AND  SOIL  DATA 

- ? - PILE-COORDINATES— A  ND-BATT-ER - 

Z  STIFFNESS  AND  FLEXIBILITY  MATRICES  FOR  THF 

STRUCTURE  AND  COORDINATES  CF  ELASTIC  CENTER 
- 4 - APPLIED -LOADS - 

5  STRUCTURE  DEFLECTIONS 

6  PILE  DEFLECTIONS  ALONG  PILE  AXIS 

- 7 - pi  l-e-forces-along-ptl-e-ax-i-s - 

P  PILE  FORCES  ALONG  STRUCTURE  AXIS 


INPUT  THE  NUMBERS  OF  THE  TABLES  FOR  WHICH  YOU  WANT  THE  OUTPUT. 
SEPARATE  THE  NUMBERS  WITH  COMMAS. 

-IM-,  2,-3;4 , - 


— IN PUT-A-F I-LE N A.  M E— F 0 R— TA B LE -8—I-N~JA— CH  A  R  AO TERS'  OR-LES-S- 
IF  YOU  WANT  TO  USE  THIS  INFORMATION  FOR  A  NEW  RUN 
HIT  A  CARRIAGE  RETURN  IF  YOU  DO  NOT  WANT  THIS  FILE. 


-Is, - 

—I N^UT-A-F IirE— NAME- FOR— OUTPUT — IN-?- CHARACTERS  OR- LES'S^ - 

HIT  A  CARRIAGE  RETURN  IF  OUTPUT  IS  TO  BE  PRINTED  ON  TERMINAL. 


I> 


EXAMPLE  PBOB  5 
—  RETAINING-WAB-t 


— NO;  -OF-PILE-ROWS— =* - 5 - B~M  A  TRT~X~ tS~~CA LCU LATED "FOR~~EACH~frOW — 


1.  TABLE  OF  PILE  AND  SOIL  DATA 
PII*  NflMBEBS 


1  5  E  =  . 15E+C7  PSI  IX  *  322.06  IN**4  IT  =  322.06  IN**4 

AREA  *  63.6  IN**2  X  =  9.00  IN  T  =  9.00  IN 

- LENGTH-*- — 30  r0-FEET - ES— -et2-.-000 - : - - - : - 

XI  *  .410?  K2  =  1.0000  K3  =  0.0000 

K4  -  0.0000  K5  =  0.0000  S6  =  0.0000 


ALLOWABLES:  COMPRESSIVE  LOAD  *  80.e00  KIPS 

TENSILE  LOAD  =  40.00e  KIPS 


TEE  B  MATRIX  FOR -PILES - -1 -THROUGH - 5-i-S- 


.778E+04  0.  0. 

0.  .265E+06  0. 

- A- - 0_ - 


2.  TAP IE  OF  PILE  COCRDINATFS  AND  5ATTER 


-  PILE '  ROW—  BATTER - U  1~C  FTt 

1  -2.40  -5.e00 

?  -2.40  -2.500 

- , - -2V40 - 0.000— 

4  VERTICAL  3.000 

«  VERTICAL  7.000 


-  3 ~ STIFFNESS-MATRIX-S-FCR-THE-STRUCTURT 


.1*?E*06  -.274S+P6  -.S22E+07 

-.P^E+PS  .  121E+C7  -.1141*08 

-.P22E*07— 114E+08  '-.324E*10 
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.1PSE-04  .480E-05  .641F-P7 

.4801-0?  .  21  P-Z-P'E  .196F-P7 

- .e<lE-f"7 - .196E-07  — -541E-09™ 


COORDI NATES  OF  ELASTIC  CENTER 

EC1  =  .f*P4  ICS  =  -.003 


-  ********  LOAD  INC  CONEITICN-1- 


•4.  MATRIX  OF  APPLIED  LOADS  0  (KIPS  S.  FEET) 
01  03  05 

- 40  .-000 - -H-e-rMie - 1737400 - 


?.  STEUCTPRE  DEFLECTION'S  <  INCHES ) 

- Dl - D? - C5 - 

-.662E-01  .e5?E-ei  .776E-03 


•>#*«■-<■«#»■  ■>**<:  <■<>#******** 


— 6-.-PI  LE-DSFL-ECTI-ONS-AL-ONO-PH:  E-AXIS— etfiCfi-ESt 
PILE  XI  X3  X5 


1 

? 

? 

4 

* 


-.102E-01  .I49E+ee  .776E-P3 

-.192E-01  .126E+00  .776E-03 

— .292E-01 - rlt?5E+00 - r776E-e3- 

- .fi62E-01  .5771-01  .776E-03 

-.«62E-01  . 2P5F-01  .776E-03 
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E44 


?.  STRUCTURE  DEFLECTIONS  (INCHES) 

.  .  j,j - D3 - 55 - 

-.?4SE+«0  .136E-C1  - ,185E-?3 

*<r *»«!***  $#$>(■<«#*«»*<■  fclfn**!**^**!***  <i <!*>*:* 


************************************************************************ 


— 8  .--PI  t£— FCRCES— A  LONO-STg«CTUR-&-*X-I-S— ffti PS-&  "Ff  £Y~) 


PILF 

FI 

F3 

F5 

1 

-16.037 

32.06° 

0.000 

2 

-16.544 

.  *  »  n  fcr>  _ 

33.329 

_ a/  cqp 

0  .000 
-  .  /*  — A  flQ- 

* 

1  n»2c - 

-2.683 

5.378 

0.000 

* 

-2.68? 

7.733 

0.000 

SUM  -55.000  113.1P0  173.400 


f.  ft  -*  #  ft  *  £  *  *  *  $  $  £  $  #  *  <•  <*  * i> *  £  *  3  :>  £  *  #  *  $ :  *  #  lie  ^  *  $  *  £  >?  *  #  $  #  Jfc  #  $  $  £  #  **  $  $  $ £  *  £  j}r  *  jfc  $  £  #  *  #: l*  * 


EXIT  • 

c>pyf 

— JCB-PP.OCFSSIKC-CCUS - 17t910- 

BY?  80/01/25.  14.46.24. 


SFLICT  DFSIHED  SIP.VICIs 
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Computations 


>  =  & 


£  aa  = 


M1  «. 

v  = 


*  5-5 
\>,  *  yi.02.' 

V- 


COMMON  ANALYTICAL  METHOD 


LONGITUDINAL  PILE  SPACING 

3'-0“  C/c 

c?Vfe.o->“4.o^  -  *2,5-  C2.S-h2.5^  >  .5* 

5 

55*+  3.0* -v  Q.52-*  C.SQ 
87*  ”75  ft2 

<2.051  +  i»3. 1  *  6*  *  27 5*1. &  "  **  -  *25o'^ 

\V3 .1  +  <L3o  d 
5  ff-J.lS’ 

22.  ta*  ±  2.62ci 

Jta  «•  3*0  ^3  *.  .5  .sW  s -2.5  -£>.5> 

“  va&  <50.48^  v>s  23.^1  *“  v*4=  \C,.e>7**  5.57k- 

113.1** ' 
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In  MJ  M  H  Ln  U>  N>  M 


Comparison  of  Results  -  Pile  Loads 


CASE  I 


COMPUTER 

ANALYTICAL 

OUTPUT  KIPS 

METHOD  KIPS 

COMMENTS 

39.1 

AO 

33. A 

3A 

27.7 

26 

Use  20T  piles 

15.3 

16 

5. A 

5.6 

CASE  II 

35.8 

AO 

Reference  Foundation  Design 

37.1 

3A 

by  Wayne  Teng.  Using  the 

38.5 

26 

analytical  method  Piles  A  6 

5. A 

16 

5  have  horizontal  load 
^1.0  :  The  Pile  batter 
for  Piles  1,  2,  3  must  be 
increased. 

7.7 

5.6 

U 


EXAMPLE  PROBLEMS  6  AND  7 


Two-dimensional  Problems 

Retaining  Wall  on  Piles  from: 
"Substructure  Analysis  and  Design" 
by  Dr.  Paul  Andersen 

Comparison  of  Computer  Output 
with  Elastic  Center  Method  and 
a  Common  Analytical  Method 


Example  Problems  6  and  7  illustrate  how  the  computer  can  be  used  to  analyze 
pile  forces  for  a  retaining  wall  founded  on  piles.  The  computer  results  are 
compared  with  the  results  obtained  by  hand  computation  methods  commonly  used 
by  civil  engineers.  The  physical  pile  layout  for  examples  6  and  7  is  shown 


Example  Problem  no.  6  compares  the  results  obtained  by  the  computer  method  with 
those  obtained  by  the  elastic  center  method  assuming  the  soil  offers  no  lateral 
support;  in  other  words,  the  subgrade  modulus  is  zero.  The  computer  results 
agree  closely  with  the  elastic  center  method  results.  A  description  of  the 
elastic  center  method  can  be  found  in  "Substructure  Analysis  and  Design"  by 
Dr.  Paul  Andersen.  This  method,  however,  is  limited  to  pile  groups  consisting 
of  hinged  piles  arranged  in  two  subgroups  whose  centerlines  intersect. 

Example  Problem  no.  7  compares  the  results  obtained  by  the  computer  method  with 
the  results  obtained  by  a  common  analytical  method.  A  description  of  the  common 
analytical  method  can  be  found  in  "Foundation  Engineering"  by  Ralph  Peck, 

Walter  E.  Hanson,  and  Thomas  H.  Thomburn,  and  in  "Foundation  Design"  by  Wayne 
Teng.  Example  Problem  no.  7  demonstrates  that  the  individual  pile  forces 
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Figure  E2.  Physical  problem  for  example  problems  6  and  7 


obtained  by  this  method  are  not  always  correct.  For  this  particular  problem, 
the  results  do  not  compare  favorably  with  the  results  obtained  by  the  computer 
method.  A  subgrade  modulus  of  312  psi  was  used  for  the  computer  method. 
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Example  Problem  6 


Two-dimensional  problem 
retaining  wall  on  piles 
no  lateral  soil  support 


Properties 

E  -  0.18  x  107  psi 

Degree  of  fixity  «  0.0 

Mod  ,  =  10  psi  0 
sub 

Pile  resistance  =  1.0 

I  -  1017.87  in.4 

X 

Participation  factor 
for  torsion  »  0.0 

I  =  1017.87  in.4 

y 

Torsion  modulus  ■  0.0 

Area  »  113.1  in.^ 

Length  =  20.0  ft 

Loading 

Case 

Q1 

(kips) 

Q3 

(kips) 

S 

(kip-ft) 

1 

-60.0 

192.0 

-1218.0 

Properties  and  loading  conditions 
for  example  problem  6 
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m 


uelco«e  to  the  hes  network 
_ sour,  access  jortis  swa_7fi. 


select  desired  service!  EKS1  (^jv) 


79/08/24.  12.44.41. 

EKS1 -175I.N0440.428— 79/08/1P.DS-0 — 02,24.12.  -79/08/2U- 
USER  NUhEERi  «S99  (g) 

-  HltllllttUtlMIIU  — - 

PASSWORD!  M  O) 

- 

IERNINAL  .  -  205.  ITT 
RECOVER/USER  10!  CHAMBERS  @) 


10  CONKENT  CN  El.S  SERVICE  CUALIir.  SEE  NSG.24. 

7HS  INrLEKENTATlON  LATE  HAS  SEEK  REVISED  ECR  LATER  IN-SEPTEH8ER. - 

USER  NOTIFICATION  Ultt  SE  TWO  WEEKS  IN  ADVANCE  Of  INPlEliENIATION. 

—  CONSOLIDATED  PLOT-LIDRART-INPLEHENTAIION  NIL4.— f£-06/24/7?  FOR-ERS1 - 

AND  EKS2.  09/09/79  ON  EKS3. 

— FOR-LAEOR  -DAT-HOUDAT  -SCHEDULE.  f'LEASE-SEE-NSS. 4. - 

- 12.45.37.COEFS_!naTNE«S/UN<EC£L8'-UF»ATED^-AUG7? - 

12.45. 37.SUBJECT  IS  CERL  FILES 

— 12.45.39*IKUCNTGRT-0F— DANS-URDATED^-AUG-79-SEE-NOTDAft/GN-CEdAT - - — 

_ C7CNEHITjSIRC.*I1_@ _ CKVggfc  /We _ jtgOMl _ 

_171  SANPLE  fP08i EH  no.  i  input  hnc.  numbers  and  c/aig 
1)2  RETAINING  NALL  ON  PILES 

— 173  2 _ _ _ 

174  5  1  1 

_  174X10 _ 

1/7  I  5  20  3 

— I>*  12 - 

I>9  2 

— 1710-2 - 

^tl  0  10  0 

_  1M2-1 - 

I>13  80  40 

-J>»4-2 - 

1715  1  2  -2 

_17i4_3-5.0 - — 

1717  2  2  5.47  7.34  13 

-  1718  -40  192.-1218  - - 

17  ©  CWO  OF  DATA  XPtPuT 

E/[Q?  -  OO-TO  TOP  OF  FltS - 

E7F  13  print  iB  lines 
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- - -  nil  ■!  rrUfartllrt  jtlfcinlliuNfi MJjwjr 


1  SAMPLE  PROBLEM  NO.  6 
-  -2  RETAINING  NALL  ON  PILES - 

3  2 

- 4  5.  LI _ 

4  1  IS 

- Z-L5-2IL3 _ 

8  12 

_ 12 _ 

10  2 

_ 1L  OJJLfi _ 

12  1 

- LL8OJ0 _ 

14  2 

_ 15X2-2 _ 

14  3  5  0 

_ IZ  2  2  5.4X134  J3 _ 

18  -40  192  -1218 

_ F>FHF  6»iT  EDIT  MOftg 

STK0N1  EDITED. 


_ DREPLACE.STR0N1  SAO£  DuOK  iTROMI 


OQLD.C  ORPS/maCtLt 
_CTCALLtCORPSjM034_J. 


program 


INPUT  DATA  FILE  NANE  IN  7  CHARACTERS  OR  LESS.  HIT  A 
_ CARKIAGEJ!EIURN_If_INFULDAIA_WILL_COHE  ERCN..IER7IINAL 


I7STR0N1  @ 


_ THIS  FROGKAM  GENERATES  .THEXOLLOHING-IRStESl. 


_ TAPLE  NO.  . 

...  -  CONTENTS 

1 

PILE  AND  SOIL  DATA 

-7 

PHF  COOPDI*ATF«  FATTfR 

3 

stiffness  and  flexibilitt  matrices  for  the 

..  ..  STRUCTURE  AND  COORDINATES  OF-  ELASTIC  CENTER.  . 

4 

APPLIED  LOADS 

5 

.  _STRuniiPF  DEFifrrin^ 

4 

PILE  DEFLECTIONS  ALONG  PILE  ARTS 

2 _ 

-  -  PILE-f  DRC£1  ALONd-PILF_AX!c 

8 

PILE  rORCES  ALONG  STRUCTURE  AXIS 

.  INPUT.  THE  NUNESES  OE. THE  TABLES  FOR  IIHICH-TSIL  NANT-THE  OUTPUT.- 
SEPARATE  THE  NUMBERS  UITH  COMMAS. 

~I>Ii2i3>4i5i4,7>J — g _ 


INPUT  A  FILENAME  FOR  TAPcE  8  IN  7  CHARACTERS  OR  LESS 

IF  YOU  NANT  TO  USE  THIS.  INFORNATION  FOR  A  HEN  RUN - 

HIT  A  CARRIAGE  RETURN  IF  TOU  DO  NOT  NANT  THIS  FILE. 
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I  I M1M4LW 


- INPUT-  A  FILE  FOR  CU7PUMN  7-  CHARACTERS-OR-ICSS* - 

HIT  A  CARRIAGE  RETURN  IF  OUTPUT  IS  TO  BE  PRINTED  ON  7EF.HINAL* 


!>  © 


_ SAMPLE  PROBLEM.  NIL— la — 

RETAINING  WALL  ON  PILES 


- NO.-OF-PJLE-RONS-- — 5 - B-NATRIX-IS-CALCULATEB-FOR-EACH-RO* 


_ iimmimmummummuuumtmuuuuuumumuttu 


1.  TABLE  OF  PILE  AND  SOIL  DATA 


PILE  NUMBERS 


.1 — 5-E-= — 18EP07-PSI — IX 
AREA  =  113.1  INI  12 
LENGTH  =  20.0  FEET 

Kl  =  .4107  K2  = 


E 


K4_? — O.COOO — K5-= 


-  —1017 . 87-IN114 — IT- *—101 7.  B7  IN1I4- 
X  =  12.00  IN  1  *  12.00  IN 

ES  =  10.000 
1.0000  K3  -  0.0000 

0.0000— K4-* — 0.0000 - 


l 


—LENGTH- DE-PILES- L  20.  OO-FEEtl-IS-INSUFtICIENT _ 

FOR  PILE  GROUP  -  1  MINIMUM  ACCEPTABLE  LENGTH  IS  38.54  FEET 

E0R-SEML-INFINITE-BEAM-0N.ELA3TIC  FOUNDATION _ 


ALLOUAELESL— COMPRESSIVE  LOAD-* - 60,000  KIPS- 

TENSILE  LOAD  *  40.000  KIPS 


THE  B  MATRIX  FGR  PILES 


1  THROUGH  5  IS 


t! ' 


.818EF03  0.  0. 

0. - .829EP04-  0.. 

0.  0.  0. 


tttiimtmimmmutijttmtutmtiiiimittitiiimimuMttim 
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2.  TASLE  Of  PILE  COORDINATES  AND  BATTER 


PILE  RON 

BATTER 

U1  (FT) 

1 - 

.2.00a 

2 

-2.00 

2.000 

3 

VERIICA1 _ 

_ 5.67Q _ 

4 

vertical 

9.340 

_ 5 _ 

-UERUCA1 _ 

- 13.000- 

_ 3 i_S  TIFENESS  . HAIR1X. S .FOR . THE. STRUCTURE. 


jJ3sEt04_uA43EL94 _ J55£L0a_ 


-.443E+04  .382EF07 

__as?E«a_^3imai_ 


.3UEFO? 

JS2EIU- 


3A  FLEXIBILITY  MATRIX  F  FOR  Th'E  STRUCTURE 


.135E-04  .454E-05  .SUE-07 

a54Er05 _ .4UE-0S _ ,  333Er07_ 

.SUE-07  .333E-07  .2R9E-0? 


COORDINATES  OF  ELASTIC  CENTER 
_ ECL_= _ .00? _ EC2_= _ .000_ 


tmmt  LOADING  CONDITION.  1  lUIIIII 


iiitmmmmmmiitmtiimmmiimiiutmmiiimmuuiu _ 


4.  MATRIX  OF  APPLIED  LOADS  0  IMPS  I  FEET) 


01  Q3  05 

-40.000  102.000  -1218.000 


I 
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f  JL£-.50RC£S-  ALONG  _SIRUCIUft£.A)(I!i_iKIES-t-EI 


-PILE _ EJ _ a _ F5- 

_1 — =29.424 - HUH - 0.004. 

2  -29.424  58,49*  0.000 

-3 - .-.,250 _ =14.250  0.000 

4  -.250  24,905  0.000 

-5 - —.250 -  43.954  0.000- 


SUM  -40.000  192.000  -1218.000 


— ittuuumuumumnmittummmmimmmmmummu 


t 

f 


EXIT.- 

_ DSTE _ k Ofi&E£ _ 

JOB  PROCESSING  CCUS  19.204 
- IYE — 29/08/24. — 12^55.43. _ 


:  select  desired  service! 

h  - - - 
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By  elosti c  center  medo d 


Consider  piles  I P  Z.  ore  pile  qroup  A 
and  p//es  3j  <f-T5  ore.  pi  I e^o  roup  6 

Then  by  the  e/orfc.  center  method 

Pa  -  P a  +.  r  M  ond  Pe,  —  Re  f  f  M 
Ha  erz  ne  sr~- 

Where  :  Rap  Re  ore  force  components  indirections  o  A 
A  f  3  pde  groups 

Ta.  C/le  ore  the  number  of  piles  m  A  f '3  pdeorojps 

r  is  the  d/stonce  from  the  efosl'C.  center  Jo 
the.  pile . 

M  >s  the  moment  of  the  oppfed  foods  about 
the  e.foz  fc  center 


Erz  =  ^6^7•)  L  -  26 .94* 
File  /  -  ZZ.oM  - . . 


Pile  Z 
Piled 

Pile  A 

Pie  5 


€2.0  M  = - 


<£) 


/g.  5 
3 


340.4)  3.67 
26-9A 


(o 

G 


ie.s  =  __ 

.  3  , 

'  L2  i_  4.  3  </■  4)  3-67 
.  3  26-9  4 


CoCo.OO 

G&.OO 
— /  3.9/ 

2  5.00 
:  63.9/ 
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Comport  son  of  resoles  obbo/ned  by  the  compofer 
me-hhod  w/th  fhe  resol/s  obbo/ned  by  bbc 
<2.1  a s he  c.en  fer  /re  fho d. 

Pile  Forces  olony  P//e  Axis 


Pile 

A Jo . 

Computer 

Outpuf 

F3 

(tips) 

E/oshc 

Center 

Fs 

(ftps) 

1 

65,147 

co.oo 

Z 

65-747 

(o<o  -OO 

3 

-14.2  SO 

-/3.9/ 

4 

24-905 

25-00 

5 

63-9  54 

63.91 
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Example  Problem  7 


TWo-dimensional  problem 
retaining  wall  on  piles 
subgrade  modulus  *  312.0  psi 
(medium  soil) _ 


Properties 

E  -  0.18  x  107  psi 

Degree  of  fixity  *  0.0 

Mod  ,  «  312.0  psi 
sub 

Pile  resistance  ■  1.0 

1  -  1017.87  in.4 

X 

Participation  factor 
for  torsion  ■  0.0 

I  -  1017.87  in.4 

y 

Torsion  modulus  “  0.0 

Area  «  113.1  in.2 

Length  *  20.0  ft 

Loading 

Case 

Q1 

(kips) 

Q3 

(kips) 

Q5 

(kip-ft) 

L 

-60.0 

192.0 

1218.0 

Properties  and  loading  conditions 
for  example  problem  7 
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io  dumi  Mia*  m  7 

20  (STAINING  WALL  ON  PILES 
302 
10  5  1  1 
SO  1  312 
60  1  5  20  3 
70  12 
80  2 
90  2 

100  0  1  0  0 
110  1 
120  80  40 
130  2 
140  1  2  -2 
150  3  5  0 

160  2  2  5.67  9.34  13 
170  -60  192  -1216 
EOI. 

07.02.03.  DARNING 


SYSTEM  DOWN  IN  5  MIN... ECS  ERRORS 
DFILE 

STROM 1  EDITED. 

DREPLACErSTROMl 
C>0LD  t CORPS/UN=CECELB 
C>CALl.C08PSrX0034 


INPUT  DATA  FILE  NAME  IN  7  CHARACTERS  OR  LESS.  HIT  A 
CARRIAGE  RETURN  IF  INPUT  DATA  NILL  COME  FROM  TERMINAL. 


DSTR0H1 


THIS  PROGRAM  GENERATES  THE  FOLLOWING  TABLES. 


TABLE  NO.  CONTENTS 

1  PILE  AND  SOIL  DATA 

2  PILE  COORDINATES  AND  BATTER 

3  STIFFNESS  AND  FLEXIBILITY  MATRICES  FOR  THE 
STRUCTURE  AND  COORDINATES  OF  ELASTIC  CENTER 

4  APPLIED  LOADS 

5  STRUCTURE  DEFLECTIONS 

4  PILE  DEFLECTIONS  ALONG  PILE  AXIS 

7  PILE  FORCES  ALONG  PILE  AXIS 

8  PILE  FORCES  ALONG  STRUCTURE  AXIS 


INPUT  THE  NUMBERS  OF  THE  TABLES  FOR  WHICH . YOU  WANT  THE  OUTPUT. 
SEPARATE  THE  NUMBERS  WITH  COMMAS. 

IMi2i3i4f5r4>7rB 


INPUT  A  FILENAME  FOR  TABLE  8  IN  7  CHARACTERS  OR  LESS 

ip  in  mr  io  m  no  amrai  m * m  mm 


i> 


EXAMPLE  PROBLEM  NO  7  * 

RETAINING  WALL  ON  PILES 

NO.  OF  PILE  ROWS  *  5  B  MATRIX  IS  CALCULATED  FOR  EACH  ROW 


nmmtmmtmmtMttmtmtummwmmutmmtmmm 


E6A 


-  -.<*«■> 


1.  TABLE  Of  PILE  AMO  SOIL  DATA 
PILE  NUMBERS 


S  E  =  .18E+07  PSI  IX  =  1017.87  1N»I4  IY  =  1017.87  IN**4 

AREA  *  113.1  INII2  X  =  12.00  IN  r  =  12.00  IN 

LENGTH  =  20.0  FEET  ES  =  312.000 

K1  =  .4107  K2  =  1.0000  K3  =  0.0000 

K4  =  0.0000  K5  -  0.0000  K6  =  0.0000 


ALLOWABLES:  COMPRESSIVE  LOAD  :  80.000  KIPS 

TENSILE  LOAD  =  40.000  KIPS 


THE  B  MATRIX  FOR  PILES  1  THROUGH  5  IS 


•108E405  0.  0. 
0.  .B29E+06  0. 
0.  0.  0. 


mtmmmmmmmtMmnnmmtmmmmmmtmtmm 


2.  TABLE  OF  PILE  COORDINATES  AND  BATTER 


PILE  RON 

BATTER 

U1  (FT) 

1 

-2.00 

2.000 

2 

-2.00 

2.000 

3 

VERTICAL 

5.670 

4 

VERTICAL 

9.340 

S 

VERTICAL 

13.000 

* 


3.  STIFFNESS  MATRIX  S  FOR  THE  STRUCTURE 


.381EW  -.A53E+OA  .157E+OG 
-.A35E404  .382E407  -.3UE409 
.157E40I  -.3UE40?  .352E411 


3A  FLEXUILITY  MATRIX  F  FOR  TME  STRUCTURE 


.7R7E-OS  .377E-05  .298E-07 
•377E-05  .273E-05  .224E-07 
.290E-07  .224E-07  .213E-09 


COORDINATES  OF  ELASTIC  CENTER 
ECI  *  .00?  EC2  *  .000 


mttm  L0MBIN6  CONDITION  i  tmtttt 


♦  m « ♦ i ♦ J7TT77 7 7 7 7 »7J7 7 7J2J777JIIJJ77J7TT±T7iJJ7 7 7 7 T 7 7 7 7 7  7  7 7 7 7 7 7 7 T T 7 7 7 7 I 


4.  MATRIX  OF  APPLIED  LOADS  0  (KIPS  1  FEET) 


01 


03 


05 


S.  STRUCTURE  DEFLECTIONS  (INCHES) 
01  03  DS 

-.103EW0  -.293E-01  -.592E-03 


6.  PILE  DEFLECTIONS  AL0N6  PILE  AXIS  (INCKS) 


PILE  XI 

X3 

X5 

1 

-.1  TIE  WO 

.684E-01 

-.592E-03 

2 

-.171E400 

.6S4E-01 

-.592E-03 

3 

-.163E+00 

.110E-01 

-.592E-03 

4 

-.103E400 

.370E-01 

-.592E-03 

5 

-.103EW0 

.630E-01 

-.592E-03 

7.  PILE  FORCES  ALONG  PILE  AXIS  (KIPS  *  FT) 

PILE  FI  H  F5  FAILURE 

SU  CO  TE 

1  *1*042  56*762  0.000 

2  -1.042  56.762  0.000 

3  -1.970  9.100  0.000 

4  -1.970  30.723  0.000 

5  -1.970  52.279  0.000 

TOTAL  NO.  FAILURES  *  0  LOAD  CASE  1 
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8.  PILE  FORCES  ALONG  STRUCTURE  AXIS  (KIPS  1  FEET) 


PILE 

FI 

F3 

F5 

1 

-27.032 

49.945 

0.000 

2 

-27.032 

49.945 

0.000 

3 

-1.978 

9.108 

0.000 

4 

-1.978 

30,723 

0.000 

5 

-1.978 

52.279 

0.000 

SUN 

-40.000 

192.000 

-1218.000 

Fy  C  ommon  Analy  hi  Co  /  /  leAhod 

Forces 

EH  -  10  k/ft 
E\l  =32  k/ft 
Mo  men  As  obouA  A 

eMa  =  io(q.s)  +  3Z(<h.o)  -  £?7  ‘V, 

Po/n.f  A  ~/~q  Resol  Aon  A 

2.17  _  8.66' 

32. 

CG  of  /?/ej  oocf  Eccentric/ Ay 

z  Cn.o )  +  /  ( 73 <0  +  /  6-67)  _  6.60  Ft  fe  PA  of  row  5 

5 

e  =  <9.66  “  <3.60  =  .  06/ 
g  onJ jSechop  MqduJij. _ 

Ed 2  *  2  6z  4  /  )  2" V  ( .  74^  +  (Z.  93)  *  +(&<sof  -  If.  44  cde-H1 

fiowj  /  ?£,  91.44  t  44/  -  20.14  pile-ft 

Ro UJ  3  7/.44-A  .74  =  /23-S7  pde-Pr 

Rolj  4  9/.  44 -A-  2.93=  3/.  2/  pile -Pi 

fio'jj  5  91*44  4  6.60  ■=  / 3.S5  pile -Pi 

Loods  or>d  MomerAs  oo  d>  Ft  Stnp> 

EH  -  to  ((a)-  &0K 
E  V  =  32  (Vj  =  /92* 

S7V-  S'/e  =  !9ZCo7>)  -  //  52  '* 
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Ver  hi  co  /  Co/npaner> ts  o  f  P/>!e  Reoc  T/on  5 
ftouj3  I  < 2.  !9Z  j  //,  52.  -  5Q.9Q>k 

5  '  20. ?4 

/?ouj  3  192  y.  II.  52  _  3a.  43  * 

5  /23-5  7 

Ro'aj  4-  I9Z  -  //.  52  •=  38.0  24 

5  31.2/ 


^  hl~C»0  kips  gH-CeO  k:ps 


P'lt  Force  Diagram  Pile  Force  D io<^rom 

Corr.puler  Method  Common  Ano/yrico!  Method 


I1 


Pile 

Forces  Along 

Pile  Axis 

Pile 

No. 

Computer 

Output 

F3  (kips) 

Comm.  Anal 
Method 

F3  (kips) 

H 

56.762 

43.56 

2 

56.762 

43.56 

3 

9.108 

38.49 

4 

30.723 

38.02 

5 

52.279 

37.57 

Comparison  of  the  results  obtained  by  the 
Computer  Method  with  the  results  obtained 
by  the  Common  Analytical  Method. 
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EXAMPLE  PROBLEMS  8  AND  9 


3-d  linens  tonal  Problems 
Aschenbrenner' s  Example 
Comparison  of  Computer  Output 
with  Aschenbrenner' s  Calculations 

Example  problems  8  and  9  illustrate  how  the  computer  can  be  used  to  analyze 
3-dimensional  pile  foundations  problems.  The  following  examples  were  taken 
from  Aschenbrenner ' s  paper,  entitled  "Three-Dimensional  Analysis  of  Pile 
Foundations",  published  in  the  ASCE  Journal  of  the  Structural  Division,  Vol. 

93,  paper  no.  5097,  Feb.  1967,  pp.  201-219.  The  computer  results  are  compared 
with  Aschenbrenner' s  calcu'  at ions.  The  physical  pile  layout  for  example  problems 
8  and  9  is  shown  in  Figure  E3. 

Example  Problem  no.  8  compares  the  results  of  the  computer  method  with  Aschen¬ 
brenner'  s  calculations  for  a  subgrade  modulus  of  35  pci,  assuming  the  subgrade 
modulus  varies  linearly  with  depth. 

Example  Problem  no.  9  compares  the  results  of  the  computer  method  with 
Aschenbrenner *8  calculations  for  a  subgrade  modulus  equal  to  zero;  in  other 
words,  assuming  the  soil  offers  no  lateral  support. 


P;/e  s 

60  P+  lony 
A~d(ol  5Qtn. 
Zx~  Q,Z10  fn * 
ly  -  &j  Z70  in* 
E-  ^ooopoo  psi 


Figure  E3.  Physical  problem  for  example  problems  8  and  9 


Example  Problem  8 


Three-dimensional  problem 
Aschenbrenner  example 
subgrade  modulus  ~  35  pci 


Properties 

E  -  3,000,000  psi 

Degree  of  fixity  “  0.0 

Modsub  “  35.0  pci 

Pile  resistance  -  1.0 

I  -  8,270  In.4 

X 

Participation  factor 
for  torsion  ■  0.0 

I  -  8,270  in.4 
y 

Area  ”  367.0  in.2 

Length  -  60.0  ft 

Torsion  modulus  “  0.0 

Loading 

Q1 

Q2 

Q3 

Q4 

Case 

(kips) 

(kips) 

(kips) 

(kip-ft) 

1 

40.0 

600.0 

0.0 

Q5 

(kip-ft) 


0.0 


Properties  and  loading  conditions 
for  example  problem  8 


}" 


10  EXAMPLE  PROBLEM  NO  Sr  ES=35 
20  ASCHENRRENNER  CHECK 
30  3 
40  9  1  1 
SO  2  35 
60  1  9  60.  2 

70  8270.  8270.  367.  16.  16. 

80  4 

90  3000000. 

100  2 

110  0.  1.  0.  0. 

120  10.  10.  10.  150.  10.  10.  150.  10. 
130  0 

140  3.  45.  5.  5.  0. 

150  3.  90.  0.  5.  0. 

160  0.  0.  5.  0.  0. 

170  0.  0.  0.  0.  0. 

190  3.  135.  -5.  5.  0. 

190  0.  0.  -5.  0.  0. 

200  3.  225.  -5.  -5.  0. 

210  3.  270.  0.  -5.  0. 

220  3.  315.  5.  -5.  0. 

230  40.  0.  600.  0.  0.  0. 

EOI. 

e^file 

8TR0M3  EDTTFO. 


C  REPLACE, STR0M3 
C>OIB,CORPS/UN=CECELB 
OCALL » CORPS  r  X0034 


INPUT  DATA  FILE  NAME  IN  7  CHARACTERS  OR  LESS,  HIT  A 
CARRIAGE  RETURN  IF  INPUT  DATA  HILL  COME  FROM  TERMINAL. 


| 


DSIMM3 


THIS  PROGRAM  GENERATES  THE  FOLLOWING  TABLES 


TABLE  HO.  CONTENTS 

1  PILE  AND  SOIL  DATA 

2  PILE  COORDINATES  AND  BATTER 

3  STIFFNESS  AND  FLEXIBILITY  MATRICES  FOR  THE 
STRUCTURE  AND  COORDINATES  OF  ELASTIC  CENTER 

4  APPLIED  LOADS 

5  STRUCTURE  DEFLECTIONS 

4  PILE  DEFLECTIONS  ALONG  PILE  AXIS 

7  PILE  FORCES  ALONG  PILE  AXIS 

6  PILE  FORCES  ALONG  STRUCTURE  AXIS 


INPUT  THE  NUMBERS  OF  THE  TABLES  FOR  WHICH  YOU  WANT  THE  OUTPUT. 
SEPARATE  THE  NUMBERS  WITH  COMMAS. 

IMi2»3r4>5»6r7*8 


INPUT  A  FILENAME  FOR  TABLE  8  IN  7  CHARACTERS  OR  LESS 
IF  YOU  WANT  TO  USE  THIS  INFORMATION  FOR  A  NEW  RUN 
HIT  A  CARRIAGE  RETURN  IF  YOU  DO  NOT  WANT  THIS  FILE. 


I> 


INPUT  A  FILE  NAME  FOR  OUTPUT  IN  7  CHARACTERS  OR  LESS. 

HIT  A  CARRIAGE  RETURN  IF  OUTPUT  IS  TO  BE  PRINTED  ON  TERMINAL. 


I> 


EXAMPLE  PROBLEM  NO  8.  ES=55 
ASCHEHBRENNER  CHECK 


NO.  OF  PILES  *  9  B  MATRIX  IS  CALCULATED  FDR  EACH  PILE 
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1.  TABLE  OF  PILE  AND  SOIL  DATA 


PILE  NUMBERS 


1  9  £  =  .30E+07  PSI  IX  -  8270.00  MU  IT  *  8270.00  IN**4 

AREA  »  307.0  IN»2  X  *  10.00  IN  Y  »  10.00  IN 

LENGTH  s  00.0  FEET  £S  =  35.000 
XI  3  .4107  K2  3  1,0000  K3  3  0.0000 

K4  3  0.0000  K'5  3  0.0000  KO  =  0.0000 


ALLOWABLES!  COMBINED  BENDING  FOR  TENSION  *  10.000  KIPS 

MOMENT  ABOUT  MINOR  AXIS  FOR  TENSION  =  10.000  KIP-FT 

MOMENT  ABOUT  MAJOR  AXIS  FOR  TENSION  3  10.000  KIP-FT 

COMBINED  BENDING  FOR  COMPRESSION  =  150.000  KIPS 

MOMENT  ABOUT  MINOR  AXIS  FOR  COMPRESSION  3  10.000  KIP-FT 

MOMENT  ABOUT  MAJOR  AXIS  FOR  COMPRESSION  *  10.000  KIP-FT 

COMPRESSIVE  LOAD  3  150.000  KIPS 

TENSILE  LOAD  3  10,000  KIPS 


THE  B  MATRIX  FOR  PILES  1  THROUGH  9  IS 


499E+05  0. 

0. 

0. 

0, 

.499E+05 

0. 

0. 

0, 

0. 

. 153E+07 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


2,  TABLE  OF  FILE  COORDINATES  AND  BATTER 


PILE  NO. 

BATTER 

ANGLE 

Ul(FT) 

U2<FT) 

U3IFT) 

1 

3.00 

45. 

5.000 

5.000 

0.000 

2 

3.00 

90. 

0.000 

5.000 

0.000 

3 

VERTICAL 

0. 

5.000 

0.000 

0.000 

4 

VERTICAL 

0. 

0.000 

0.000 

0.000 

S 

3.00 

133. 

-5.000 

5.000 

0.000 

L 

VERTICAL 

0. 

-5.000 

0.000 

0.000 

7 

3.00 

225. 

-5.000 

-5.000 

0.000 

8 

3.00 

270. 

0.000 

-5.000 

0.000 

9 

3.00 

31 S. 

5.000 

-5.000 

0.000 

3.  STIFFNESS  MATRIX  S  FOR  THE  STRUCTURE 


.745E+04  -.212E-03 

.749E-03 

-.922E-01  -.753E+08 

-.B94E-07 

-.212E-03  .104EM7 

0. 

.1296+09  -.5406-01 

-.191E-01 

.7A9E-03  0. 

. 129E+08 

-.477E-04  .477E-M 

.1466+00 

-.922E-01  A&IM 

-.477E-06 

.2986+21  0. 

-.439E401 

-.753EF08  -.540E-01 

.477E-06 

0.  .3096+11 

.387E-04 

-.894E-07  -.191E-01 

. 14AE+00 

-.4396+01  .3B7E-04 

•215E+10 

3A  FLEXIBILITY  MATRIX 

F  FOR  THE 

STRUCTURE 

. 4 786-05 

-.1966-15 

-.1066-15 

.6356-17 

.4346-08 

-.409E-23 

-.1946-15 

•204E-05 

-.328E-21 

-.8866-08 

.3126-17 

.1846-18 

-.1066-15 

-.3286-21 

.7776-07 

.2656-23  -.2596-18 

-.527E-I7 

.6356-17 

-.8866-08 

.2656-23 

.7176-10  -.9706-26 

.675E-19 

•434E-08 

•312E-17 

-.2596-18 

-.1796-25 

.4306-10 

-.592E-24 

-.+096-23 

.1846-18 

-.5276-17 

.6756-19  -.5926-24 

.444E-09 

tmm*  LOADING  CONDITION  1  ttmttt 
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4.  MATRIX  OF  APPLIED  LOADS  0  (KIPS  t  FEET) 


01  02  03 
40.000  0.000  000.000 


04  OS  00 

0.000  0.000  0.000 


5.  STRUCTURE  DEFLECTION  (IJKHES) 

01  D2 

D3  D4 

D5 

DO 

.713E-01  -.784E-11 

.400E-01  .254E-12 

.174E-03 

-.310E-11 

ttmtmummmmtummmmmMttmmnmnmmmtm 


PILE 

XI 

X2 

a 

X4 

X5 

U 

1 

.304E-01 

-.504E-01 

.S03E-01 

.117E-03 

. 123E-03 

.389E-04 

2  - 

.147E-01 

-.713E-01 

.442E-01 

•105E-03 

-.310E-12 

.549E-04 

3 

•713E-01 

-.198E-09 

.302E-01 

.254E-12 

.174E-03 

-.310E-11 

4 

.713E-01 

-.784E-11 

.400E-01 

.254E-I2 

.174E-03 

-.310E-11 

5  - 

.059E-01 

-.504E-01 

.382E-01 

.117E-03 

-.123E-03 

.389E-04 

0 

.713E-01 

.iea-w 

•570E-01 

.254E-12 

•174E-03 

— .310E-11 

7  - 

■•059E-01 

.504E-01 

.382E-01 

-,117EH)3 

-.123E-03 

-.38PE-04 

8  - 

-.I47E-01 

•713E-01 

.442E-01 

-.105E-O3 

.441E-12 

-.549E-04 

9 

.304E-01 

.504E-01 

.503E-01 

-.U7E-03 

.123E-03 

-.389E-04 

7.  PILE  FORCES  ALONG  PILE  AXIS  (KIPS  1  FEET) 


PILE 

FI 

F2 

F3 

F4 

F5 

F6 

CBFTR 

I 

1.814 

-2.514 

76.858 

0.000 

0.000 

0.000 

.51 

2 

-.735 

-3.555 

67.606 

0.000 

0.000 

0.000 

.45 

3 

3.555 

-.000 

55.321 

0.000 

0.000 

0.000 

.37 

4 

3.555 

-.000 

71.263 

0.000 

0.000 

0.000 

.48 

5 

-3.284 

-2.514 

58.354 

0.000 

0.000 

0.000 

.39 

A 

3.555 

.000 

87.205 

0.000 

0.000 

0.000 

.58 

7 

-3.284 

2.514 

58.354 

0.000 

0.000 

0.000 

.39 

8 

-.735 

3.555 

67.606 

0.000 

0.000 

0.000 

.45 

9 

1.814 

2.514 

76.858 

0.000 

0.000 

0.000 

.51 

FAILURE 
CB  BU  CO  TE 


TOTAL  NO.  FAILURES  =  0  LOAD  CASE  1 


8.  PILE  FORCES  ALONG  STRUCTURE  AXIS  (KIPS  t  FEET) 


PILE 

FI 

F2 

F3 

F4 

F5 

FA 

1 

20.180 

16.625 

72.340 

0.000 

0.000 

0.000 

2 

3.555 

20.682 

64.369 

0.000 

0.000 

0.000 

3 

3.555 

-.000 

55.321 

0.000 

0.000 

0.000 

4 

3.555 

-.000 

71.263 

0.000 

0.000 

0.000 

5 

-9.068 

12.623 

56.398 

0.000 

0.000 

0.000 

6 

3.555 

.000 

87.205 

0.000 

0.000 

0.000 

7 

-9.068 

-12.623 

56.398 

0.000 

0,000 

0.000 

8 

3.555 

-20.682 

64.369 

0.000 

0.000 

0.000 

9 

20.180 

-16.625 

72.340 

0.000 

0.000 

0.000 

SUN 

40.000 

-.000 

600.000 

0,000 

-.000 

.000 

Pile  Forces  Along  Pile  Axis 


Computer 
Output 
F3  (kips) 


76.858 


67.606 


55.321 


71.263 


58.354 


87.205 


Aschen- 
brenner 
F3  (kips) 


79.7 


67.5 


49.6 


71.3 


55.3 


93.0 


58.354 


55.3 


67.606 


76.858 


67.5 


79.7 


Comparison  of  the  results  obtained  by  the 
Computer  Method  with  the  results  obtained 
by  Aschenbrenner. 
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Example  Problem  9 


Three-dimensional  problem 
Aschenbrenner  example 
no  lateral  soil  pressure 


Properties 

E  -  3,000,000  psi 

Degree  of  fixity  =  0.0 

Mod  ,  =  .  lpci  «  0 
sub 

Pile  resistance  =1.0 

I  =  8,270  in.4 

Participation  factor 

for  torsion  =  0.0 

I  =  8,270  in.4 

Torsion  modulus  =  0.0 

y 

Area  »  367.0  in.^ 

Length  «  60.0  ft 

Loading 

Case 

Q2 

(kips) 

m 

*4 

(kip-ft) 

Q5 

(kip-ft) 

Q6 

(kip-ft) 

1 

40.0 

0.0 

600.0 

0.0 

0.0 

0.0 

Properties  and  loading  conditions 
for  example  problem  9 


E82 


10  EXAMPLE  PROBLEM  NO  0.  £S=0 
29  ASCHEN8RENNER  CHECK 

39  3 

40  9  1  1 
50  2  .1 

*0  1  9  40.  2 

70  8270.  8270.  347.  14.  14. 

80  4 

«0  3000000. 

100  2 

110  0.  1.  0.  0. 

120  10.  10.  10.  150.  10.  10.  150.  10. 
130  0 

140  3.  45.  5.  5.  0. 

150  3.  90.  0.  S.  0, 

140  0.  0.  5.  0.  0. 

170  0.  0.  0.  0.  0. 

180  3.  135.  -5.  5.  0. 

190  0.  0.  -5.  0.  0. 

200  3.  225.  -5.  -5.  0. 

210  3.  270.  0.  -5.  0. 

220  3.  315.  5.  -5.  0. 

230  40,  0.  400.  0.  0.  0. 

EOI. 


.  .ai1' 

DREFLACE.STR0H3 
D0LD.C0RPS/UN-CECEL8 
C>CALL » CORPS  *  X0034 


INPUT  DATA  FILE  NAME  IN  7  CHARACTERS  OR  LESS.  HIT  A 
CARRIAGE  RETURN  IF  INPUT  DATA  HILL  COME  FROM  TERMINAL. 


I7STR0N3 
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THIS  PROGRAM  GENERATES  THE  FOLLOWING  TABLES 


TABLE  NO.  CONTENTS 

1  PILE  AND  SOIL  BATA 

2  PILE  COORDINATES  AND  BATTER 

3  STIFFNESS  AND  FLEXIBILITY  MATRICES  FOR  THE 
STRUCTURE  AND  COORDINATES  OF  ELASTIC  CENTER 

A  APPLIED  LOADS 

5  STRUCTURE  DEFLECTIONS 

b  FILE  DEFLECTIONS  ALONG  PILE  AXIS 

7  PILE  FORCES  ALONG  PILE  AXIS 

8  PILE  FORCES  ALONG  STRUCTURE  AXIS 


INPUT  THE  NUMBERS  OF  THE  TABLES  FOR  WHICH  YOU  WANT  THE  OUTPUT. 
SEPARATE  THE  NUMBERS  WITH  COMMAS. 

I>l»2f3f4f5f4f7f8 


INPUT  A  FILENAME  FOR  TABLE  8  IN  7  CHARACTERS  OR  LESS 
IF  rOU  RANT  TO  USE  THIS  INFORMATION  FOR  A  NEW  RUN 
HIT  A  CARRIAGE  RETURN  IF  YOU  DO  NOT  WANT  THIS  FILE. 


I> 


INPUT  A  FILE  NAME  FOR  OUTPUT  IN  7  CHARACTERS  OR  LESS. 

HIT  A  CARRIAGE  RETURN  IF  OUTPUT  IS  TO  BE  PRINTED  ON  TERMINAL. 


I> 


EXAMPLE  PROBLEM  NO  9<  ES=0 
ASCHENBRENNER  CHECK 


HO,  OF  PILES  *  9  B  MATRIX  IS  CAECULATED  FOR  EACH  PILE 


PILE  NUMBERS 


1  9  E  =  . 30E+07  PSI  IX  *  8270,00  INM  IY  =  8270.00  INM 


AREA  = 

347.0 

m»2 

X  = 

14.00  IN  Y  •  16.00  IN 

LENGTH 

=  60, 

0  FEET 

ES  = 

.100 

K1  = 

.4107 

K2  = 

1.0000 

K3  = 

0.0000 

K4  = 

0.0000 

K5  = 

0.0000 

K4  = 

0.0000 

allowables:  combined  bending  for  TENSION  =  10.000  KIPS 

MOMENT  ABOUT  MINOR  AXIS  FOR  TENSION  =  10.000  KIP-FT 

MOMENT  ABOUT  MAJOR  AXIS  FOR  TENSION  =  10.000  KIP-FT 

COMBINED  BENDING  FOR  COMPRESSION  =  150,000  KIPS 

MOMENT  ABOUT  MINOR  AXIS  FOR  COMPRESSION  =  10.000  KIP-FT 

MOMENT  ABOUT  MAJOR  AXIS  FOR  COMPRESSION  =  10.000  KIP-FT 

COMPRESSIVE  LOAD  =  150.000  KIPS 

TENSILE  LOAD  =  10.000  KIPS 


THE  B  MATRIX  FOR  PILES  1  THROUGH  9  IS 


.148E+04 

0. 

0. 

0. 

.14BEE04 

0. 

0. 

0. 

•153EE07 

0. 

• 

0. 

0. 

0. 

• 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

0. 


tmmmtmmmmmmmmtmmmmtmummmmmm 


2.  TABLE  OF  PILE  COORDINATES  AND  BATTER 


PILE  NO. 

BATTER 

ANGLE 

Ul(FT) 

U2(FT) 

U3(FT) 

1 

3.00 

45. 

5.000 

5.000 

0.000 

2 

3.00 

90. 

0.000 

5.000 

0.000 

3 

VERTICAL 

0. 

5.000 

0,000 

0,000 

4 

VERTICAL 

0. 

0.000 

0.000 

0.000 

5 

3.00 

135. 

-5.000 

5.000 

0.000 

t 

VERTICAL 

0. 

-5.000 

0.000 

0.000 

7 

3.00 

225. 

-5.000 

-5.000 

0.000 

8 

3.00 

270. 

0.000 

-5.000 

0.000 

9 

3.00 

315. 

5.000 

-5.000 

0.000 

tttt1ttttt1ttUni*t1t*ittttt*tttl****f****+***************»******tt*** 


3.  STIFFNESS  MATRIX  S  FOR  THE  STRUCTURE 


•319ET0& 

-.219E-03 

.794E-03 

-.952E-OI 

-.778E408 

-.154E-06 


-.219E-03 

,624Em 

-.373E-08 

.133EI09 

-.558E-01 

-.197E-01 


. 794E-03 
-.373E-08 
•128EF08 
0. 

-.477&.06 

.151ET00 


-.952E-01 

,133EW 

0. 

.297E411 

0. 

-.453E40I 


-.778E+08 
-.558E-0J 
-,477E -06 
0, 

.308E+11 

.396E-04 


-.154E-04 

• 15IE+00 
-.453E+01 
.396E-04 
.641E+0B 


3A  FLEXIBILITY  MATRIX  F  FOR  THE  STRUCTURE 


.815E-05 

-.149E-13 

-.504E-15 

.102E-15 

.206E-07 

.485E-20 


-.149E-13 
.318E-04 
.922E-20 
-.142E-06 
• 14SE-16 
-.253E-15 


-.504E-15 

.922E-20 

.778E-07 

-.412E-22 

-.127E-17 

-.183E-J5 


.102E-15 

-.142E-04 

-.412E-22 

.667E-09 

-.324E-23 

.351E-17 


•206E-07 

.148E-14 

-.127E-17 

-.345E-23 

.843E-10 

-.272E-23 


.485E-20 

-.253E-15 

-.183E-15 

.351E-17 

-.272E-23 

.I56E-07 


UUtM  LOADING  CONDITION  1  HUM* 


4.  MATRIX  OF  APPLIED  LOADS  0  (KIPS  1  FEET) 

01  02  03  04  05  Q6 

40.000  0.000  600.000  0.000  0.000  0.000 


nmmmmntmmmmtmmtmmmmtmmmmmmtm 


5.  STRUCTURE  DEFLECTIONS  (INCHES) 

D1  02  03  04  05  06 

.326E+00  -.678E-09  .467E-01  .407E-11  .823E-03  -.110E-09 


mmmtmmummmummmummmmmmmmmmtt 


6.  PILE  DEFLECTIONS  ALONG  PILE  AXIS  (INCHES) 


PILE 

XI 

X2 

X3 

X4 

X5 

X6 

1 

.220E400 

-.2J1E400 

.704E-01 

• 552E-03 

.582E-03 

. 184E-03 

2 

-.148E-01 

-•326E+00 

•443E-01 

.780E-03 

-.437E-11 

.260E-03 

3 

•326E400 

-.728E-08 

-.265E-02 

•407E-11 

.823E-03 

-.I10E-09 

4 

.326E400 

-.678E-09 

.467E-01 

.407E-11 

.823E-03 

-.110E-09 

5 

-.249E+00 

-.231E400 

•I82E-01 

•552E-03 

-.582E-03 

. 184E-03 

6 

. 326E+00 

.593E-08 

•961E-01 

.407E-11 

.823E-03 

-.II0E-09 

7 

-.249E400 

.231E400 

. 182E-01 

-.552E-03 

-.5B2E-03 

-.184E-03 

B 

-.148E-01 

•326E400 

.443E-01 

-.780E-03 

.496E-11 

-.260E-03 

9 

•220E400 

•231E400 

. 704E-OI 

-.552E-03 

.502E-O3 

-.184E-03 

tttmummmnmmmnmmmmmmumumtmnumm 
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7.  PILE  FORCES  ALONG  PILE  AXIS  (KIPS  l  FEET) 


PILE 

FI 

F2 

F3 

F4 

F5 

F8 

CBFTR 

1 

,226 

-.342 

107.650 

0.000 

0.000 

0,000 

.72 

2 

-.022 

-.484 

87.758 

0.000 

0.000 

0.000 

.45 

3 

.484 

-.000 

-4.048 

0.000 

0.000 

0.000 

.40 

4 

.484 

-.000 

71.424 

0.000 

0.000 

0.000 

.48 

5 

-.370 

-.342 

27.888 

0.000 

0.000 

0.000 

.19 

8 

.484 

.000 

148.893 

0.000 

0.000 

0.000 

.98 

7 

-.370 

.342 

27.888 

0.000 

0.000 

0.000 

.19 

8 

-.022 

.484 

87.758 

0.000 

0.000 

0.000 

.45 

P 

.326 

.342 

107.850 

0.000 

0,000 

0.000 

.72 

TOTAL  NO.  FAILURES  *  1  LOAD  CASE  1 


tmtmmmmmmmtmtmmmmmmmmmtttmtmtm 


8.  PILE  FORCES 

ALONG  STRUCTURE  AXIS 

(KIPS  1 

FEET) 

PILE 

FI 

F2 

F3 

F4 

F5 

F  8 

1 

24.532 

24.048 

102.023 

0.000 

0.000 

0.000 

2 

.484 

21.408 

84,288 

0.000 

0,000 

0,000 

3 

.484 

-.000 

-4.048 

0.000 

0.000 

0.000 

4 

.484 

-.000 

71.424 

0.000 

0,000 

0.000 

5 

-5.741 

8.225 

28.553 

0.000 

0.000 

0.000 

8 

.484 

.000 

148.893 

0.000 

0.000 

0.000 

7 

-5.741 

-8.225 

28.553 

0.000 

0,000 

0.000 

8 

.484 

-21.408 

84,288 

0.000 

0,000 

0.000 

9 

24.532 

-24.048 

102.023 

0.000 

0,000 

0.000 

SUN 

40.000 

.000 

800.000 

.000 

-.000 

,000 
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Pile  Forces  Along  Pile  Axis 


Pile 

No. 

Computer 

Output 

F3  (kips) 

Aschen¬ 

brenner 

F3  (kips) 

1 

107.650 

112.5 

2 

67.758 

67.5 

3 

-4.046 

r-13.1 

4 

71.424 

71.3 

5 

27.866 

22.5 

6 

146.893 

155.7 

7 

27.866 

22.5 

8 

67.758 

67.5 

9 

107.650 

112.5 

Comparison  of  the  results  obtained  by  the 
Computer  Method  with  the  results  obtained 
by  Aschenbrenner . 


ana**  MMumfre 


In  accordance  with  latter  froa  DAEN-RDC,  DABN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Card*  for 
laboratory  Technical  Publications,  a  facsinila  catalog 
card  in  Library  o t  Congress  NOC  fornat  is  reproduced 
below. 


CASK  Teak  Group  on  Pile  Foundations. 

Basic  pile  group  behavior  /  by  the  CASK  Task  droop  on 
File  Foundations.  Vicksburg,  Hiss.  :  0.  B.  Waterways 
Experiaeirt  Station  ;  Springfield,  Va.  ;  available  frxm 
Rational  Technical  Infomation  Service,  1900. 

22,  [ill]  p.  :  ill.  ;  27  cn.  (Technical  report  - 
U.  S.  Amy  bgineer  Waterways  fspertaant  Station  ;  K-fiO-5) 
Prepared  for  Office,  Chief  of  bgineers,  D.  8.  Arsqr, 
Washington,  D.  C. 

References:  p.  21-82. 

1.  Cbaputer  progress.  2.  Oonputerlsed  siaulation.  3.  Design 
criteria.  1.  Pile  foundation  design.  5.  Pile  foundations. 

I.  United  States.  Amy.  Corps  of  fciglneers.  II.  Series: 
United  8tates.  Waterways  Xxperlaant  Station,  Vicksburg, 

Hiss.  Technical  report  ;  1-80-5. 
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UM  of  Computer  Programs  for  Computer-Aided 
Structural  Engineering 


Instruction  Report  0-79-2  User’s  Guide:  Computer  Program  WHti  interactive 
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